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ABSTRACT 


The plane strain problem of adhesively bonded structures which 
consist of two different - isotropic adherends is considered. By 
expressing the x-y components of the displacements in terms of 
Fourier integrals and using the corresponding boundary and 
continuity conditions, the system of integral equations for the 
general problem is obtained. Then , these integral equations are 
solved numerically by applying Gauss-Cbebyshev integration scheme. 

The shear and the normal stresses in the adhesive are 
calculated for various geometries and material properties for a 
stiffened plate under uniaxial tension <r x . Also the numerical 
results involving the stress intensity factors and the strain energy 
release rate are presented. The closed-form expressions for the 
Fredholm kernels are provided, so that the solution for an arbitrary 
gecmetry and material properties can easily be obtained. 

The numerical solution of the integral equations indicates that 
as (h^/a), (hjj/a) and (hg/a) decrease the convergence becomes slower 
and hence computations become costlier. For the general geometry, 
the contribution of the normal stress is quite significant. For the 
symmetric geometries, however, the dominant stress is the shear 
stress. More specifically, the normal stress vanishes if the 
adherends also happen to be of the same material and the same 
thickness. 





1. INTRODUCTION 


In order to optimize performance and fuel consumption, 
aerospace and marine industries have been turning to the use of 
advanced (fiber- reinforced organic) composites, more and more in 
commercial aircraft, military aircraft and marine systems. These 
materials offer very good strength- to- weight and stiffness-to-weight 
ratios. However, one major drawback is the strength and fatique 
penalty introduced by mechanical fasteners at joints. So more 
sophisticated joining methods are required. 

Adhesive bonding, on the other hand, provides a desirable 
alternative to mechanical fastening because of; 

*1. Load being carried over a larger area, thus reducing the 
stress concentration, 

2. Higher joint efficiency (relative strength-to-weight of 
the joint region), 

3. No decrease in strength due to fastener holes, 

4. Less expensive and simpler fabrication techniques, and 

5. Lower maintenance costs. 

However, adhesive bonding has its awn disadvantages. The load 
is not carried uniformly over the entire bond area, but instead is 
confined to a small region along the bond edge. Though not as high 
as the stresses at a rivet, this highly stressed region, can lead to 
failure. 


2 



The past forty years have witnessed the expenditure of 
considerable analytic effort in an attempt to describe stress-strain 
distributions in composite structures formed by. the adhesive bonding 
of materials. The efforts of Goland and Reissner [1], have been 
extended by the computerised and experimental analyses of numerous 
investigators. 

To gain some insight and to provide criteria for further 
development of bonding materials and bonding techniques, assumptions 
have been introduced which are justified only by the analytic tools 
available to the investigator. Goland and Reissner [1], for-example, 
restrict themselves to adherends of the same material having 
identical length and thickness, with no stress variation within the 
adhesive film. With progress in analytic techniques, each succeeding 
investigator has been able to relax the number of assumptions 
previously required to obtain a solution. 

However, because of the nonhomogeneous nature and of the 
geometrical complexity of the medium, even for the linearly elastic 
materials the exact analytical treatment of the problem regarding 
the stress analysis of the structure is, in general, hopelessly 
complicated. The existing analytical studies are, therefore, based 
on certain simplyfying assumptions with regard to the modeling of 
the adhesive and the adherends. The adherends are usually modeled as 
an isotropic or orthotropic membrane [2], a plate [1,3, 4, 5] or an 
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elastic continuum [5,6]. The adhesive on the other hand, is usually 
treated as a shear spring [2,6], a tension shear-spring [1,7], or is 
neglected [8]. In this report the adhesively bonded joint problem is 
considered by assuming both the adhesive and the adherends as 
elastic layers. 

In order to design adhesively bonded structures with high 
degree of reliability, one needs to recognize that their failure 
mode is characterized by flow growth and progressive crack 
propogation. 

The energy balance criterian for fracture, based on works of 
Griffith [9] and Irwin [10] is adopted. It supposes that fracture 
occurs when sufficient energy is released from the stress field to 
generate new fracture surfaces at the instant of crack pro paction. 
This strain energy release rate provides a measure of the energy 
required to extend a crack over a unit area, and is termed the 
fracture energy. In this report, the fracture energy of an adhesive 
layer will be determined, since this property has been widely 
recognized as the appropriate criterion for adhesive failure as in 
[5,11,12,13,14]. 
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2. FORMULATION OF THE PROBLEM 


2. 1 Equilibrium Equations 


The problem considered is a stiffened plate shown in Figure 1 , 
under the following assumptions; 

- The medium is can posed of homogeneous, isotropic, elastic 
layers with different mechanical properties, 

- The problem is one of plane strain, that is, the bonded 
joint is very "wide", 

- The only external load acting on the medium is the 
uniaxial tension, cr^ x = 9 Q away from the reinforcement 
region. 

In the plane theory of elasticity the equations of equilibrium 
in terms of displacements for the isotropic materials can be 
expressed as; 


+ \t) + nV 2 u + X = 0 , 

c)x 


(?v + |#) + v + Y = 0 , 


( la, b) 


where u,v are the x, y- components of the displacement vector, X-Y are 
the x-y components of the body force vector, 
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and P are the shear modulus and the Poisson's ratio, 

respectively. 

For each of the layers shown in Figure 1 , and for no body 
forces, the equations (1a,b) read as, 

+ h ± \7 2 u = 0 , 

+ H i y 2 v = 0 , i=1 ,. . ,4 (2a, b) 

2.2 Solutions u ± , v ± , <T* y 

As it is seen from Figure 1 , the medium possesses a geanetric 


c) e 

( + Hi )—— 

Ox 


C)e 

( + H ± ) 

c)y 


symmetry with respect to x=0 plane, so the problem is solved for 
x >, 0 . Also note that, the x-y components of the displacements may 



be expressed as Fourier integrals, since the displacements as well 
as their derivatives decrease sufficiently rapidly as |x(— >oo , so 
that the requirement of absolute integrability is satisfied. 

Therefore, assuning the x-y components of the displacements in 
the i'th layer in the form, 


u ± (x,y) 



(a,y) Sin(Ctx) da 


v-j/x.y) 




(a,y) Cos(ax) da 


o 


(3a, b) 


and using the field equations (2a, b), one obtains, 


2 r°° 

u i (x,y) = — l[ (A^ + A^y) e“ ay + 

•n 


(a^ + A il} y) e ay 3 Sin(ax) da 



o 



Cos(ax) da 


(4a, b) 
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where A i j ' s are functions of CC which will be determined from the 
continuity and the boundary conditions. Afterwards, the stresses are 
evaluated by Hooke's Law, and expressed as, 


.00 


yy 

2M, v J 


-[ctCAfi + A^y) + 2( 1-«? ± ) A^ J e’ uy 


^-fltfA^ + A^y) + 2(1-1^) A^ J e ay j Cos(CCx) 


da , 


.00 


^yx = 

2M ± J ir 


[-[ a( A^ + A^y) + (1-2H i ) A ^2 3 e _CCy 


[ aCA^ + A^y) - (1-21^) A ±i} ] e^y j Sin(ax) da 


(5a, b) 


2.3 The Boundary and the Continuity Conditions 

On the boundaries y=h^ , y=-h !| , the medium possesses the 

following homogeneous boundary conditions; 

^^y(x,h., ) = 0 , 0 < x ( 00 (6a) 
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''iy^’V = 0 

9 

0 < x <oo 

(6b) 

ffjy(x.-h^) = 0 

9 

0 ^ x <co 

(6c) 

*Jy(X’-V = 0 

9 

0 4 x <oo . 

(6d) 


The continuity conditions require that on the interfaces the 
stress and the displacement vectors in the adjacent layers be equal, 
that is, 


tf^y(x,h 2 ) - fl^x.l^) _ o , 

llyU.V - -T^yCx.l^) = 0 , 

u 1 (x,h 2 ) - u 2 (x,h 2 ) = 0 , 

(x, h 2 ) - v 2 (x,h 2 ) = 0 , 

«-^ y (x,-h 3 ) - ff3 y ( Xf _ h3 ) = o f 

^^y(x,— h^) _ ^^y(x,“ h^) =0 , 

u^Cxj-h^) - u^Cxj-hg) = 0 , 

V|j(x,-h 3 ) - v 3 (x,-h 3 ) = 0 , 


0 4 x <co 
0 4 x < oo 
0 4 x <00 
0 ^ x < oo 
0 4 x <00 
0 ^ x < oo 
0 4 x <oo 
O s < X <00 


(6e) 

(6f) 

( 6g) 

(6h) 
(6i) 
(6 j) 
(6k) 
( 61 ) 
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The above conditions (6a-l) provide 12 linear homogeneous 
algebraic equations in terms of 16 unknowns. So 4 more equations are 
needed. Those are obtained from the surface which has the crack, 
that is, at y s 0 , 


ff^y(X,0) ~ = 0 

, 0 ^ x < oo 

(7a) 

*2y(x,0) - *5cy(x, 0) = 0 

, 0 ^ X < 00 

(7b) 

^yy^x, 0) = tf^yCx.O) = g(x) 

, x£L 

(8a) 

T^y(X,0) = (X, 0) = f(x) 

, x£l . 

(8b) 


L is the part of the x-axis without the crack and f(x), g(x) 
are respectively, shear and normal stresses at the very same region. 
The mixed boundary conditions at y=0 , and the process of 

superposition as shown in Figure 2, give rise to the integral 
equations for the problem. Those are, 


lin y-*o 

a 

Six 

[ u 2 (x,y) 

- u 3 (x,y) ] 

= * , 

0 ^ X < 00 


a 

Ox 

[ v 2 (x,y) 

- v 3 (x,y) ] 

= 0 . 

0 N < x < oo (9a,b) 
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Note that the integral equations have been expressed in terms 
of the first derivatives of the displacement differences with 
respect to x. Also note, X appearing in equation (9a), has the 
following values depending on the geometry of the medium, 

for plane stress 
for plane strain . 

2.4 Application of the Boundary and the Continuity Conditions 
In equations (8a) and (8b), it has been assumed that 

0 , x > a 

f-l (x) , x ^ a 

0 , x > a 

f 2 (x) , x 4 a ( 1 Oa, b) 

where "a" is the bond length as shown in Figure 1. Also note ‘that 
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equations (5a, b) at y=0 gives, 





^xyU.O) = 

+ [«a 23 - 



+ (1-2F 2 )A 22 ] 
Sin(Otx) da , 


0 4 x <oo 


0 4 x <oo . (11a, b) 


Ihe above equations with the conditions stated in (10a, b), read as, 


f 2 (x) 

2M 2 


2 r°° 

~ JV 


F 2 (a) Cos(ax) da 


0 4 x ^ a 


1 2 r 00 

f^x) = — I F.j(a) Sin(ax) da , 0 4 x 4 a (12a, b) 

2M 2 « q J 


wi th, 



F 2 (cc) = - CCA 21 - (1-2? 2 )A 22 + aA 23 - (1-2H 2 )A 21j . (12c, d) 

Conditona (6), (7) and (12 c,d) provide 16 algebraic equations 
to be solved for the same number of unknowns in terms of and F 2 . 
Writing these equations in matrix form, 




A 11 

1 

1 


' 0 ' 

» 

i 

Coefficient 


« 

1 

a 

| 


1 

1 

1 

i 

Matrix 


A ij 

= 

0 

1 

f( a , h ±t V ± ) 


i 

i 

i 

i 


1 

1 

F^Ct) 



A 2}4 


1 

0 

u - 


(16x16) (16x1) (16x1) 

and multiplying each side with the inverse of the coefficient 
matrix, gives A^'s. 




1 | 

< 1 


-i 

0 

i 

A ij 

= 

Coef. ' 

,Q ij 

1 

s u 


i 

0 



, I 

Matrix 1 

• 1 

• ; 



I 

F^CC) 



i ! 



f 2 («) 
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where Q^'s and S ij' 3 (i»J = 1 - ,4) are the 15'th and 16* th coloumns 
of the inverse coefficient matrix, respectively. After performing 
the matrix multiplication, 






it may easily be shown that A.y's can be solved in the following 
general form, 


A ij (CC) = Q i j(0)F 1 + s i;j (a)F 2 . (16) 


Uiese when substituted in (13) give two systems of equations to 
solve for Q^j ' s and S^j's as follows, 



0 


6 

0 

1 
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which leads to; 







'o ' 

1 

Coef. 





1 

1 

1 

Matrix 



Q ij 

s 

1 

1 

0 






1 

. 

- 


- 


0 

- 

- 


- 


’ 0 ’ 

1 

Coef. 



S U 


1 

\ 

1 

1 

Matrix 




0 






0 






1 


( 18 ) 


(19) 


(18) and (19) are solved first numerically for every desired 
value of CC . It is definitely much easier and less time consuming 
process rather than trying to solve them analytically. However, it 
has its own shortcomings. It has to be kept in mind that certain 
combination of these and S^j's (envelope functions) will 

actually be integrated from zero to infinity at every x and t (see 
Section 2.5-6). This requires (18) and (19) to be solved at 
sufficiently many CC's. This is obviously a very costly and time 
consuming job. Especially for the thinner geometries, where the 
convergence of the envelope functions is very slow, (18) and (19) has 
to be solved at even more CC's, in order to achieve certain 
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significant figure accuracy in the evaluation of the infinite 
integrals. 

So as a result, (18) and (19) are required to be solved 
analytically. After rather lengthy manipulations, the closed-form 
expressions for Q^'s and S^'s are found. 

Note that, it is most convenient to retain Q 22 and Q 2 4 as the 
final two unknowns in equations (6a, b, e, f, g,h) and (12c,d). 
Equations (12c, d) actually provide two equations for Q^’s and 
Sf-j’s, namely; 

— C(A 2 ^ — 2( 1-? 2 )A 22 — CCA 2 2 + 2( 1-H 2 )A 2 jj = F^(CC) 

- CCA 21 - (1-2J» 2 )A 22 + «A 23 - ( 1-2*^ 2 )A 24 = F 2 (0) 

by substituding (16) into above equations we find, 

— OC Q 2 i — 2 ( 1 -P 2 ) Q 22 “ + 2 ( ‘•“•'2)024 = ^ 

- «Q 21 - (1-2H 2 )Q 22 + 0CQ 23 - (1-2P 2 )Q 24 = 0 (20a, b) 

and 

- aS 21 - 2(1 -i> 2 )S 22 - «S 23 + 2(1-P 2 )S 22j = 0 

— 0CS 2 ^ "* ^ ^”^^2)^22 "* ^^23 ” ( 1~2P 2 )S22 j = 1 . (21a, b) 
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Solving (6a, b,e, f, g, h) , (20a, b) one gets; 


G 3 G 5 “ G 2 G 6 

°2 G 4 - G 5 G i 

g 1 g 6 ” G 3 G i} 

GgGij - 


Q 21 = ~ZZ [ _1 " K 2 Q 22 + Q 2il] 
2CC J 

Q 23 = — -[l ~ Q 22 + K 2 Q 24 ] 


( 22 ) 


(23) 


(24) 


(25) 


Q 14 = [- 


-2CC£(c-1) 


1+K. 


■3 Q 21 + [■ 


-2a£(c-l)(h 2 + k 2 / 2a) 

1-tK, 


3 Q. 


22 


+ C 


0 + K, 


1 + K. 


3 0,4 


( 26 ) 


Q i2 = C- 


c + K, 


1 + K. 


3 Q 22 + [■ 


-2a(i-c) 
fid + k 1 ) 


3 Q 


23 


-2a(i-c)(h 2 - k 2 /2a) 

+ [ ] Q, 

fid* «,) 


24 


(27) 


M Kj 

G 13 = [— — ] Qi 2 + [ - h 2 + — ] Q- 4 

2 a 2 a 
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( 28 ) 


+ [ — J Q 22 + $23 + ^ ^2 ** — ' ^ 

2 fiC 2 a 


K 1 -*1 

Q = [ - h 2 ] Q 12 + [ - 

20 2a* 


] Q 


14 


k 2 k 2 

+ Q 2 -j + C ^2 + — 3 ®22 + ^ -I ^24 * 

2a 2a* 


Similarly, retaining Q32 and Q^ij in equations (6e,d,i, 
is the easiest way to obtain, 


*32 


G 14 G 9 " G 1 2 G 1 1 
G 1 0 G 12 " G 13 G 9 


^34 


G 13 G 11 ~ G 14 G 10 
G 10 G 12 “ G 13 G 9 


^31 


= “ 77 b + K 3 Q 32 " Q 34 ] 

2CC 


Q 33 = “ b “ Q 32 + k 3 Q 34 ] 


K, + d -2a* (1-d) 

Ql}2 = [__£ ] q 32 + [ ] Q 


1 + K, 


1 + K, 


33 


+ [ 


-2a* (d-i)(h 3 + k 3 /2a) 


■] Q 


1 + K, 


34 


( 29 ) 


,k,l), (7a, b) 


( 30 ) 


( 31 ) 


( 32 ) 


( 33 ) 


( 34 ) 
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Q 44 = [ 


-2«(d-1) 

*\l + Kj.) 


3 Q 3 1 + [ 


d + K- 


1 + K, 


■3 Q. 


34 


-2a(d-1)(-h 3 + K 3 /20) 

+ [ | ] 

* (1 + K k ) 


( 35 ) 


Q41 = ( b 3 “ — 1 3 Qjip + C 

2 a 




2 a 


3 Q. 


44 


1C jEj ^ 

+ Q 3 . + [-h, + -'■ i — 1 3 Qop + ( *~*~ 3 Qoji 

J 0 2a 3 2a J 


(36) 


K4 

Qii 3 = [ ;3 Qmo + [bo + 3 Qji]i 

5 2a£ 5 2a 


-K, k 3 

+ [ _ 3 Qqp + Q 33 + — ■ 3 Qqij , 

2a£ 5 J 2a J 


(37) 


G, , (1=1 ,16) , *, »* are defined in Appendix II. 

Similarly using equations (6a, b, c, d, e, f ,g,h, i, j,k, 1) 
(21a, b), one obtains, 


P 3 P 5 " P 2 P 6 
P 2 P 4 - P 5 P 1 


(38) 


and 
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s 24 = 


P 1 P 6 " P 3 P 4 
P 2 P 4 - P 5 P 1 


(39) 


1 

3 21 = ~ C " 1 " K 2 S 22 + ^24 3 
2CC 


(40) 


23 


1 

[-1 — Sop + ^?^24 ^ 

2a 


(41) 


*14 


-2a£ (c-1) -2a£(c-1)(hp + Kp /2a ) 

[ ] s 21 + [ ] S 22 


1 + K 


1 


1 + K 


1 


+ c- 


C + K, 


1 + K 1 


•] 


(42) 


S 12 = C 


K 2 + c 

1 + K, 


-2a(i-c) 

'] Sop + ( ————— ] Sop 

*(1 + K,) 


-2a(l-c)(hp - Kp /2a ) 

+ [ ] s^ 

*(1 + K,) 


(43) 


JBk- k. 

S-o = [ S«p + (“hp + — — ] S 

- 2 a 20 


14 


-£l < 2 K 2 

+ C— -3 s 22 - s 2 + [h 2 - — ] s 

2a 2a 


24 


(44) 


K 

S n = C-(h 2 + — ) ] S 12 + [ — ] S u 

2a 2a* 
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p i3 p ii ~ p i4 p io 
P 10 P 12 " P 13 p 9 

p 14 p 9 " P 12 P 1 1 
P 10 P 12 “ P 13 P 9 

1 

- ■ C 1 + K3S32 — 3 


_ [ -1 _ s 32 + k 3 s 34 ] 

20c 


Kn + (1 

l— ] s, 2 + [ 

1 + K„ 


-2a % ( 1-d) 


1 + Ki 


3 S, 


33 


+ [ 


-2a£ (d-1)(h 3 + k 3 /2a) 

1 + k., 


3 s. 


34 


-2a(d-1) 

fi ( 1 + K), ) 


3 s 3 i + [ 


d + K. 


1 + K. 


-3 S. 


34 


-2a(d-1)(-h, + k, /2a) 

+ C * 3 S^ 2 

* (1 + KJ 


( 45 ) 


(46) 


(47) 


(48) 


(49) 


(50) 


( 51 ) 



"41 


[h 3 - 


2a 


J + C” 




] S 


2a 


44 


JB*k, 


+ S31 + [-hg + ] S3? + [ 

2 a J 2a 


•] S 34 


(52) 


^43 = ^ S42 + [ hg + ' ] Sjjjj 


- k 3 * 3 

+ [ "3 S^p + ^35 + (“^3 “ J S?j, 

2 a* J J:5 J 2 a 


(53) 


P i (a) , (i=1, ,16) , Jt , * , c, d are defined in Appendix III . 


2.5_ Integral Equations 

Integral equations for the problem are derived from (9a, b), 
that is, 


o r 1 

lim y _* 0 — |^u 2 (x,y) - u 3 (x,y)J = , 0 ^ x < 00 


lim y _ 0 — [v 2 (x,y) - v 3 (x,y)j =0 , 0 ^ x <co 
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where, one might recall 


X = 


for plane stress 


* = 




for plane strain 


Substituding (4a) into above equations, and then replacing the 
corresponding Ajj's with (16), the following equations are obtained, 


.03 


lim 


y— >o 


'IT 


O 


[(Q 21 (a) + Q 22 (cc)y) (Q 23 (a) + Q 24 (a)y) e ay ] 


- [(Q 31 («) + Q 32 (cc)y) e _cty + (Q 33 (a) + Q 34 (a)y) e ay ] 


. F 1 (OC ) Cos(CCx) da 


.03 


llm y^o- 

ir 


o 


[(S 21 (a) + S 22 (a)y) e“ ay + (S 23 (a) + s^tajy) e ay ] 


- [(S 31 (a) + s 32 (a)y) e- ay + (S^a) + S 34 (a)y) e ay ] 


. F 2 (a) Cos(ax) da = 7 \ 


(54) 
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.CO 


-a 


[( Q 21< a > + Q2 2 (K 2 /« + y)) e _CCy 

+ (-Q 2 3 (ct) + /flt “ y 0 e<Xy 

(Qj^ec) + Q 32 («)(K 3 /cc + y)) e‘ ay 


(-Q 33 (a) + Q 3 J| («)(k 3 /a - y)) e ay 


. F.,(a) Sin(OCx) da 


+ lim 


y—o 


m 


.00 


-a 


(s 21 (a) + S2 2 (a)(K 2 /a + y)) e- ay 


(-s 23 (a) + s 24 (a)(K 2 /a - y)) e ay 
(^^a) + s 32 (a)(K 3 /a + y)) e“ ay 
^-s 33 (a) + s 34 (a)(K 3 /a - y)) e* y 


• F 2 (a) Sin(ax) da = 0 
(55) 


F.j , F 2 must be replaced by Fourier inversion of (12a, b) 
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F^CC) = 


J 

o 


.00 

f,(t) 

2»o 


Sin(CCt) dt 


f 2 «x) = 


.00 


j 

o 


f 2 (t) 

2» 2 


Cos(at) dt 


(56 a, b) 


Noting that, 


fyt) = f 2 (t) = 0 , t£L* 


y = 0 (57) 


where L' is the part of the x-axis containing the crack. Therefore, 
(55a, b) with (57) would read as follows , 


F 1 (a) 


F 2 (a) 


a 



f^t) 

2» 2 


a 



f 2 (t) 
2 » 2 


Sin(cct) dt 


Cos(ctt) dt . 


(58a, b) 


Replacing F^a) and F 2 (0t) in ( 5-4) , (55) with (57a, b) give two 
equations of the form ; 




a 

a 


Ih, (x,y,t) f,(t) dt = 2JJ 5 X 
j-l J J * 


o 
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11m — * 

y— *o 

o 


a 

Xh 2 J (x f y,t) fj(t) dfc = 0 . 


( 59 a, b) 


Here, 


■/ 


co 


h'j-jCxjyjt) = I RE 1 (a,y) Cos(CCx) Sln(Ott) da 


J 


CO 


h 12 (x,y,t) = J TE 1 (a,y) Cos(ax) Cos(at) da 


■/ 


oo 


b 2 -|(x,y,t) = | RE 2 (a,y) Sin(ax) Sin(at) da 
o 

r°° 

h 22 (x,y,t) = I TE 2 (a,y) Sin(ax) Cos(at) da 


( 60 a) 


( 60 b) 


( 60 c)' 


( 60 d) 


where, 


RE 1 ( a , y) = a 


[ ( Q 2 1 ( a ) + Q 22 («)y) e“ ay + (Q 23 («) + Q 24 (<X )y) e ay ] 
- [(Q 31 (a) + Q 32 (a)y) e _ccy + (Q 33 (a) + Q 3 ]J (CC)y) e ay ] 


TE 1 (a ,y) = a 


[(S 21 (a) + S 22 (a)y) e“ oy + (S 23 (a) + ^(ajy) e ay ] 
- [(s^ta) + s 32 (a)y) e- ay + (S 33 (a) + S 34 (a)y) e ay ] , 


RE 2 ( a , y) = -a 


Q 21 (a) + q, 2 («)(k 2 /a + y) j e" ay 
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+ 



+ q 24 (cc)(k 2 /« 

+ Q 32 («)(»«2 /Ct 

+ Q 31 < («)(k 3 /a 


y) ] e«y J 
y) ] e“®y 
y) ] e ay 


TE2(CC, y) 



+ S2 2 («)(k 2 /« + y) ] 
+ s^(«)(k 2 /a - y) ] 
+ s 32 (a)(K 3 /a + y) j 
+ s 34 (a)(K 3 /a - y)] 


e-«y 



e-«y 

e«yl . 


(6la,b,c,d) 


For the physical problem under consideration, the displacement 
differences on the bond surface are known and the stresses, f t ) 
are unknown, which may be determined frctn the integral equations 
given by (59a, b). After determination of fj_(t), all the desired 
quantities, like the stress intensity factor, the strain energy 
release rate and stresses can easily be evaluated. 

Here it should be clearly noted that in deriving (59a, b), the 
derivatives of the displacements, rather than displacements 
themselves, have been used. Also note that the integral equations 
should be solved under the following equilibrium conditions 
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f 

a 

I 


a 

f t (t) dt 
a 

f 2 (t) dt 


= 0 


= 0 . 


(62a, b) 


The infinite integrals giving and h 22 can be expressed as a 
sum of two integrals, the integrands of which, respectively, are of 

0(e" ay ) and 0(e _2Gch 2, e“ 2flCh 3) for OC ►oo. The first leads to a 

Cauchy Kernel and the second to a Fredholm Kernel. On the other 
hand, the integrands of the infinite integrals giving h^ 2 and h 21 

are of 0(e -2ah 2, e -2flt ^3) for a ^oo , hence leading to Fredholm 

Kernels only (see Section 2.6} . 


2.6 Cauchy and Fredholm Kernels 


After performing asymptotic expansion as a — *-co» G i *s and P^s 
(Appendix II, Appendix III) are found to be ; 



28 



G 6 = a 2 

G 16 = 0 

, P 6 = -a 2 

G ? = 0 


P 7 = o 

G 8 = a 3 a 2 » 


P 8 = ~ a 3^2 » 

G 9 = -2a , 


P 9 = -2ab 4 , 

G 10 = b 2 » 


P 10 = b 2 


where (a^, i=1,..,4), (b i , 1=1,. .,4) are given in Appendix I. Using 
them in (22) ,(23) ,(24) ,(25), (30) ,(31), 132) ,(33) and 
(38) ,(39) ,(40) , (41) , (46) ,(47) ,(48) ,(49) gives the following results 
as <X — fc-oo, 


Q 21 (a) 

= -(1 + k 2 ) / 2 a 

of 

N 

ft 

II 

O 

• 

O 

'Q 22 (cO 

= 1.0 

Q 32 (a) = 0.0 

q 23 («) 

= 0.0 

Q 33 (a) = (1 + k 3 ) / 2a 

Q 2i| («) 

= 0.0 

Q 34 (a) = 1.0 , 



(63a, b, c, d, e, f,g, h) 

and - 



S 2 i(a) 

= (-1 + k 2 ) / 2 a 

S^ta) = 0.0 

S 22 ( ® ) 

= - 1.0 

^gfa) = 0.0 

s 23 («) 

= 0.0 

S 33 («) = (-1 + k 3 ) / 2a 
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Sg^CC) = 0.0 


S 34 (CC) = 1.0 


(64a,b,C,d,e,l*,g,h) 


Substituting them into equations (64a,b, c,d), and noting K3 =k 2 , 
give , 


RE1(cc— *oo, y + ) = REloo 


-(1 + K ? ) 


2 fit 


+ y - 


(1 + k 3 ) 

2a 


+ y 


3 -ay 


REloo = [ -(1 + K 2 ) + 2ay ] e _cty , 


(65a) 


TE1 ( a — ~oo,y + ) = TEloo 


r -(1 + 


= a 


- y - 


2« 


-1 + K. 


2a 


+ y 


3 -ay 


TEloo = 0.0 


(65b) 


RE2(a— » 00, y + ) = RE2co 


= -a 


-(1 + Kp) 


2a 


— + y + 

a 


(1 + k 3 ) 

2a 


■- y 


,-ay 


RE2co= 0.0 


(65c) 
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TE2(<X 


,y + ) = TE2 


r “1 + Kr 


-a 


-1 + K- 


- y + 


2 a 


2 a 


- y 


»~ay 


TE2oo= [ (1 + Kp) + 2a y ] e’ ay 


( 65 d) 


It is obvious that if the infinite parts are substracted from 
the integrands, the integrals will be uniformly convergent and give 
bounded kernels, that is, Fredholm kernels, so limy_^ 0 , can be put 
under the integration sign. We then obtain, 


R1(a ,0) = lim y _^ 0 [ RE1(a,y) - RE1 m (a,y)] 

= cc [ Q 21 (a) + Q 23 (a) - Q 31 (a) - Q 33 (a) ] + (i + k 2 ) , 

T1(a,0) = lim y _ 0 [lE1(a,y) - TE1«(a,y)] 

= a [ S 21 (a) + s 23 (a) - s^a) - s 33 (a) 3 , (66a, b) 

R2(a,0) = lim y __ 0 [ RE2(a ,y) - RE2®(a,y)] 

k 2 k 2 

= -a [ Q 21 (a) + — Q 22 (a) - Q 23 (a) + — Q 2 4 (ct) 
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K K 

-Q 31 («) Q 32 (o) + Q3 3 ( <X ) ^Q 34 («) ] , 

cc cc 


T2(O,0) = lim y _ 0 [ TE2(Ct , y) - TE2(a,y)] 


K- K 2 

= “C [ ^21 (<X) + ■ " $ 22 ^ ^ " ^23^^^ + ~ ""^24^^^ 


K K 

- S^tcc) - 3 S 3 2 (0t) + S 33<°> -^(a) ] 

CC cc 


- ( 1 + Kg) . (66c r d) 


Integral equations (59a, b) with (66a,b,c,d) would then become ; 


lim. 



co 

RE1o=(a,y) Cos(CCx) Sin(CCt) dCC dt 


2 r a 00 

+ __ / r- (t ) / hi (a) Cos(ax) Sin(at) da dt 

v J o 



00 

T1(a) Cos(ax) Cos(at) da dt = 2w 2 * 


(67a) 
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iiaiy 


o* 



oo 

TE2m(Ct ,y) Sin(Ctx) Cos(at) da dt 




oo 

T2(a) Sin(ax) Cos(at) da dt 



00 

R 2(a) Sin(ax) Sin(at) da dt = 0. 0 


(67b) 


In (67a, b) the kernels in the first terms are Cauchy type 
whereas the remaining kernels are bounded in the closed interval 
0 ^ (x,t) ^ a . Next we change the integration limits from (zero-to- 
a), to (-a to +a). Frcm the physics of the problem, we note that 
^(Xjy) is odd, <T(x,y) is even, therefore, 


f^t) = -f^-t) 


-a < t < +a 


f 2 (t) = f 2 (-t) 


-a < t < +a . (68a, b) 


and the integral equations (67a,b) become 
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lim y— o* — 


f^t) 


+ 

-a 


RE1 (a, y) Cos(Ox) Sin(«t) da dt 


j f\,(t) K^Cx.t) dt + J f 2 (t) K 12 (x,t) dt = 2M 2 * 


’ J 

-a 

j 

o 

a 

f^t) 

K 11 


.a 


Ilia 


y-H»o + 


IT 


/ 


00 


- a 
a 


f 2 (t) J TE2 (a,y) Sin(ax) Cos(at) da dt 
o 


+ J f^t) JC 21 (x,t) dt +J f 2 (t) K 22 (x,t) dt = 0. 0 . 

(69a, b) 


-a 


_a 


where the fredholm kernels are given by, 


: n (x,t) = — J 

<tr J 


co 


Ki i (x, t) =— | R1(a) Cos(ax) Sin(at) da , 

m 


. 00 


K 12 (x,t) = — Ti(a) cos(ax) cos(at) da , 


.00 


K 21 (x,t) = — R2(a) Sin(ax) Sin(at) da , 


ir 
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1 


00 


K 22 (x,t) = 


T2(CC) Sin(OCx) Cos(«t) dCt 


(70a, b, c, d) 


The first parts of (69a, b) are reduced to Cauchy integrals as 
follows ; 


lin y-c 


-(1 + K 2 ) 




+ a 

r ,00 




f,(t) e" ay Sina(t-x) da dt 


-(1 + K,) 


+a 

r 


lim 


y-*.o + 


j 

-a 


(t-x) 

f 1 (t) 2 1 

(t-x) 2 + y 2 


dt 


+a 


-(1 + k 2 ) [ f,(t) dt 

ir J (t-x) 

-a 


(71a) 


and 


+a 


.00 


( i + k ) r — 

limy — f 2 (t) e“ ay (-Sina(t-x)) da dt 


-a 


+a 

+ k 2 ) f (t-x) 

lim y-o + f 2 (t) 5 o dt 

* y J (t-x) 2 + y 2 

-o 
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-(1 + K ,) 


+ a 
r 


j 

-a 


f 2 (t) dt 
(t-x) 


(71b) 


2.7 Normalization of the Integral Equations 

To solve the integral equations, it is convenient to define the 
following dimensionless variables and functions ; 


r = x/a 

9 

-a 4 x 4 +a 

s = t/a 

9 

-a ^ t ^ +a 

. 6 = oca 

9 

0 < a <00 

g^s) = 

f 1 (s. a) 

/ *0 

g 2 (s) = 

f 2 (s. a) 

/ *0 • 


(72) 


With (72), (69) and (62) can be written as, 


1 


+1 

r 


g^s) ds 
(s-r) 


/ 


+ I g-j(s) K 11 (r,s) ds 


+ ! 


+ I 82 ( 3 ) K^ 2 (rjS) ds = P^j 


(73a) 
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(73b) 


with, 


and 


f 


g^(s) da = 0. 0 , 


1 g 2 (s) da 

« J (s-r) 


+ J g 2 ( 3 ) K 22 (r, a) da 

-I 
+ 1 

+ j g-|(s) K 21 (r,a) da = 


0. 0 


(73c) 


with, 


+1 


/ 

-1 


g 2 (a ) da = 0. 0 


(73d) 


Here, 


.00 


K-jl(r,3) = - 


K 12 (r,s) = - 


K 21 (r, a) = - 


K 22 (r , s) — ■■ 


ir(1 + 

k 2 ) 

1 


ir(l + 

k 2 ) 

1 


u(l + 

k 2 ) 

1 


ir( 1 + 

k 2 ) 


R1(S/a) Coa(8r) Sin(Sa) d8 , 


.00 


T1(S/a) Coa(Sr) Cos(8a) dS , 


.00 


R2(8/a) Sin(8r) Sin(8a) d& , 


oo 


T2(S/a) Sin(Sr) Coa(8a) d& , 

(74a,b,c,d) 
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and, 


2H 2 (1 - »jj) 
e 4 (i + k 2 ) 


2.8 Evaluation of the Infinite Integrals in the Fredholm Kernels 


Ihe integrands of the infinite integrals (74b, d) go to infinity 
as 8-~0, that is, they have a pole at 8=o . Hence, those integrals, 
if treated separately, will be divergent and their evaluation 
requires special care. 

As 8 approaches to zero, Cos(8s) goes to unity, therefore the 
integrands become independent of "s" thus, because of the single 
valuedness condition, (73d), the coefficient of the unbounded 
integrals would vanish. Equations (74b, d) could then be replaced by 
> 


oo 

If 

K 1P (r,s) = T1(8/a) Cos(8r) (Cos(8s) - 1.0) d8 , 

*( 1 + Ko) J 


K 22 ( r 


.=> = 

1T(1 + K 2 ) J 


CD 


T2(8/a) Sin(8r) (Cos(8s) - 1.0) d8 


(75a, b) 


38 



3. SOLUTION OF THE SINGULAR INTEGRAL EQUATIONS 


3.1 Solution of the Integral Equations 


The solution of the singular integral equation 


+ 1 


r 



+ k(x,t) ♦(t) dt = g(x) , 

t-x J 


-1 < x < +1 


(76) 


subject to the single valuedness condition 

+1 


/ 


♦ (t) dt = 0 




(77) 


is given in [15]. The method has been summarized in Appendix IV . 

However, the singular integral equations (73a, c), that appear 
in this report, have two unknown functions g^ and gg. So the 
solution method that has been described in Appendix IV, should be 
modified accordingly. Also note that, since g i (s) has a power 
singularity 1/2 at the end points, the solution will be sought in 
the form ; 


gi (s) = (1 - s 2 )" 1/2 ♦ i (s) 


(78) 
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where ^(s) is a function defined in the interval -1 ^ s < +1 and 
the indices of the singular integral equations are +1. 

Following the procedure described in [13]» the integral 
equations may be expressed as 


IT 1 N-l 

* I f 

— k 11 (r k ,s 1 ) ^(s^ +2, k n (r k' s i) 

N-1 2 


1 # 

+ — k n (r k ,s N ) 
2 


IT 

M-1 


1 M-l 

f T-* # 

— k 12 ( r k»s'i) + L k 12^ r k» s ’i^ ^2^ 3 ' i^ 

2 i-2. 


+ - k 12 (r k' s ’M ) ♦2 (S 'M ) 


= P 


1 » 


k=1 , ( N— 1 ) (79a) 


with 


1 n-1 1 

— *l( 3 l) + I ^(Si) +“ ♦^Sn) = 0 > (79b) 

2 i=a, 2 


i*0 



and 


*TT 

N-1 


1 S-l 

| T-« | 

— k 21 (rj.a,) ♦ l(3l) +2- k 21^ r j> s i^ ♦i(Si) 


i-a 


+ — ^(r^Sjj) 

2 


M-1 


1 M-l 

§ y’* t 

— k 22^ r j» 3 'l^ $2^ 3 ' 1 ^ + f- k 22^ r j» s 'i^ ♦2^ s 'i' 


v-a. 


+ 7 


= 0 


j=l (M-1) 


with 


1 M-l 1 

— 02(3*1) +X ^ 2 ^ S ' +— ^2 ^ s 'm^ = 0 

2 Ua. 2 


where 


k* 1 (r, s) =— + Ki 1 (r, s) , 

ir s-r 


a 

ki 2 (r,s) = Ki 2 (r,3) 


ill 


(80a) 


(80b) 




(r,s) 


= K 21 (r,s) 


k 22 (r,s) 


1 1 

— + K 22 (r,s) • 

v s-r 


( 8la,b,c,d) 


3he kernels K,,(r,s), K, 2 (r,s), K 21 (r,s), K 22 (r,s) are defined 
by (74a, 75a, 74c and 75b), respectively, and, 


3 i = Cos 


1-1 


N-1 


it 


i=1 ,. . . . ,N 


r k = Cos 


2k-1 


2N-2 


k=1 ,....,( N-1 ) 


s' ^ = Cos 


1-1 


M-1 


m 


, i- 1 i • ••• ,M 


r j = 003 


2J-1 


2M-2 


, j=1 ,. ... , ( M— 1 ) . (82a,b,c,d) 


Equations (79a, b) and (80a, b) would give N+M unknowns 


♦ ,(31) , 1=1 . ,N 


and 


♦3(3'^ , 1=1, ,M . 
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3.2 The Stress Intensity Factors 


In Section (2.J») , it has been assumed that 


^(x.O) = f i (x ) 

-a ^ x 4 +a 

<Tyy(X, 0) = ^(X) . 

-a < x < +a 


They have been non-dimensionalized in Section (2.7), by 
defining the following variables and functions, 


r = x/a , -1 ^ x ^ +1 

(r) = f.j (r. a) / « , 0 , 
g 2 (r) = f 2 (r. a) / «‘ 0 , 


where g^s themselves are defined in the following form (see section 
3.1), for the solution of the integral equations; 


S 1 (r) = (1-r 2 )“ 1/2 ^(r) , 
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<p 2 (r) • 


g 2 (r) = (1-r 2 r U2 


In plane strain crack problems, the symmetric and anti- 
symmetric components of the stress intensity factors may be defined 
as, 


k 1 = lim x _ a -^2(a-x) C^(x f 0) , (83) 


k 2 = llmjj^ -\J 2(a-x) <T X y(x, 0 ) . (84) 


Substituting the above definitions in (83) and (84), the 
constants k^, k 2 may be related to the functions $ 2 as follows: 


k-| = lim r _.j yiTn-r) f 2 (a. r) , 

= lim r _ i _ i y2a(1-r) g 2 (r) , 


. $ 2 (r) 


= lim r ^i V2a «T ( 


* 2 (r) 


V* 


+r 
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k 1 = ff oV* ♦2( r=1 ) 


(85a) 


similarly, 


kg = llm r _^ 1 ^2a(1-r) f.,(a. r) 

= lim r _ i>1 *y/2a(1-r) S Q g.,(r) 

I »i^? 

* 1J “r*1 V 2a <'- r ) ®o 


11m 


r-^1 


-fil 


♦ ■j(r) 


Vi 


+r 


k 2 = *oV* ♦ 1 (r=1) . 


(85b) 


So the normalized stress intensity factors are 


*i ■ —h= = ♦2 (r=,) 

» 0 v a 


k 2 = — = ♦ 1 (r=1 ) ( 86a, b) 


where, ♦ 1 (r=1), ♦ 2 (r=1) are 
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$1 (s^, i*1 *— *■ S., = 1 ) | 

♦ 2 (s'i, i=1 — ► s» -,=1 ) , (87) 

can be obtained from the solution of integral equations (79) and 
(80) . 


3.3 The Strain Energy Release Rate 


The strain energy release rate is defined as, 


G = 


Kf + k| 


( 88 ) 


where, 


E = E , for plane stress 

E = E / ( 1 - P 2 ) , for plane strain 

K =Vir k 
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So for our problem, which is a plane strain case, we have 


where k^'s 


Note that 



are defined by (85a, b), one gets 


G = 


1 + K 2 

1 fla *2 ($2(1) + +§(1)) 

8M 2 


( 89 ) 


♦ ^'s are given by equations (87). 
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4. RESULTS AND DISCUSSION 


The system of integral equations, with the additional 
conditions (73 a,b,c,d) is solved by the technique described in 
Appendix IV and Section (3.1). The main problem encountered is the 
evaluation of the infinite integrals rather than the solution of the 
integral equations. Unless these integrals, which are to be used 
later to set up the algebraic system, calculated with sufficiently 
high accuracy, the solution of the integral equations will not 
obviously produce dependable results beyond certain number of 
significant digits. So the numerical integration scheme has a vital 
importance. 

Before going into any integration technique, one should first 
be able to define the envelope functions, R£1(<X), TEI(Ct), RE2(Ct), 
TE2(CC), (61 a, b,c,d). Previously, these functions (as mentioned in 
Section 2.4) were unavailable in closed form. All of their desired 
properties, including the asymptotic behaviours for <X approaching 
infinity or zero, were accurately determined by plotting and using 
certain tricks, such as multiplying with the powers of alpha in 
order to find out the singular behaviour around zero. Since the 
analytical expressions for the envelope functions were not 
available, they had to be calculated numerically, by solving sixteen 
algebraic equations each time, say "n" times, for every CC^ 
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The number "n", which has been mentioned in the previous 
paragraph, is actually the number of points needed for the numerical 
integration scheme to give sufficient accuracy. So, obviously "n" 
depends on the convergence of the integrands in (74 a,c) and (75 
a, b). More specifically, slower the convergence, larger will be the 
"n", or vice versa. A closer look, however, reveals that (74c) and 
(75a) do not in themselves, have a convergence problem. The 
convergence of (74a) and (75b), on the other hand, is greatly 
affected by the adhesive thickness. As it can be seen from Table 1 ,2 
, convergence is very slow for thin geometries. So, as a result, we 
end up with a very costly procedure, trying to solve sixteen 
equations "n" times. Because of this major drawback, equation 
systems (18) and (19) are solved analytically in order to obtain 
closed-fora solutions for the envelope functions. 

After having defined the envelope functions, one can then 
proceed with the selection of the appropriate integration scheme for 
each of the infinite integrals (74 a, c) and (75 a,b). As mentioned 
in the previous paragraph, the integrands of the equations (74c) and 
(75a) do not possess a convergence problem. So, they are integrated 
from zero to infinity by using Laguerre polynomials in one step. 
However, the other two, that is, (74a) and (75b) needed special 
treatment, because of the behaviours of their integrands. 



After a closer look, one realizes that the envelope functions 
of these integrands (74a, 75b) can be handled more easily and 

accurately, if they are examined in three intervals, instead of just 
one going from zero to infinity. First, from zero to A , where the 
functions are very steep and relatively large in magnitude. Then 
from A to B , where the functions are smoothly decreasing in 
magnitude, and finally from B to infinity where the envelope 

functions have been replaced by their asymptotic expressions for 
large alpha for computational reasons. The constants A and B , 
depend on the geometry and the material properties of the medium, 
however, in general A lies between 1 and 5 » and B is around 500 . 
So the equations (74a), (75b) are integrated in three steps using 
Filon’s integration scheme, first from zero to A, then A to B, and 
finally from B to infinity. 

Finally, the system of integral equations (73 a,c) is solved 
using Gauss-Chebyshev integration scheme as discussed m Appendix 
IV. It should be pointed out that, in order to build up the 
algebraic system for the solution, the kernels, K.,.j(r,s), K 12 (r,s), 
K 2 -|(r,s), K 22 (r,s), (74 a, c and 75 a,b), that is, the infinite 

integrals, have to be evaluated (N x (N-1)), (M x (N-1)), (N x 

(M-1)) and (M x (M-1)) times at corresponding "r" and "s n , 

respectively. N and M are the number of points for which the unknown 
functions ♦ 1 (s J _), 1=1,. ..,N , ♦2^ s j^> j=1,. ..,M , are evaluated (see 
Section 3*1). However, taking into consideration the symmetry and 
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the anti- symmetry of the kernels and the unknown functions and 
♦ 2 f these numbers can be reduced to 1/4 * th of their original values 
which means a great deal of saving in computational time and money. 
This concludes almost all of the numerical considerations regarding 
the solution of the problem. The results obtained are discussed in 
the following paragraphs. 

The material constants used in the calculations, unless 
otherwise is specified, are as follows, 

Upper and Lower Adherends : Aluminum 


E 1 = E 4 = 10.5 x 10 6 psi , 

Vi = = 3.9474 x 10 6 psi , 

= P4 = 0.33 , 

Adhesive : Epoxy 

E 2 = E 3 = 0.28 x 10 6 psi , 
p 2 = M3 = 1.0 x 10 5 psi , 

9 2 - l>3 * 0.40 . 

The results for the specimen with equal thickness adherends 
which has the same material properties have been tabulated in Tables 
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3-10. Corresponding graphs are Figures 3-6 for Tables 4,5,7 and 8 
respectively. Since the only applied load is uniaxial tension, the 
specimen is free to bend. Consequently, in the case of identical 
adherends the normal stress in the adhesive is found to be zero, 
which agrees with [4]. Also, it was observed that the adhesive shear 
stress, the corresponding stress intensity factor and the strain 
energy release rate increase as the adhesive thickness decreases. 
Table 5 shows the effect of adherend ticknesses on the adhesive 
shear stress. 

A special geanetry is studied in Tables 9 and -10, for 

com pari sion with [4]. If the stresses are to bb calculated at 

specific distances away from the right end, rather then at specific 
values of the non-dimensional variable (x/a), the similarity will 
become apparent. So we may conclude that the stresses are 

independent of bond length, hence, the strain energy release rate 
turn out to be constant (Figure 7). The similar result is found in 
[4] by using tne plate theory. 

Tables 11-20 and Figures 8-10 give the results for the specimen 
having similar adherends with different thicknesses. In Table 14, 

upper plate is less stiff than the lower plate, while in Table 11 
the relative stiffness is reversed. This is accomplished simply by 
varying the adherend thicknesses. The peak normal stress changes 
from tension in Table 14 to compression in Table 11, while its 
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magnitude remains almost the same. The shear stress is compressive 
in both cases as expected. Tables 13 and 15 shows the effect of 
(h^/a) and (hjj/a) ratios on “the adhesive stresses, respectively. - 
Again for h^ > h^ (Table 13) there is compressive normal stress, for 
h.| < hjj (Table 15) there is tensile normal stress. Also it has been 
observed that the adhesive stresses, the stress intensity factors 
and the strain energy release rate increase as the adherend 
thickness increases. The same behaviour can also be seen in Tables 5 
and 8 . This trend has been noted in [3], [ 6 ]. 

Tables 21-28 give the results for the specimen having 
dissimilar adherends. The adhesive shear stress increases as the 
shear modulus of the upper plate, Hi increases relative to the shear 
modulus of the lower plate h^ » no matter what their thickness ratio 
is. However, the peak normal stress is compressive and increases 
with increases with increasing Hi for hi > h^ (Table 21), tensile 
and decreases with increasing Hi for hi < h^ (Table 22), and is 
tensile for Hi < Hij» zero for Hi = compressive for Hi > H 4 for 
equal thickness adherends (Tables 23,24). Same trend is observed for 
the corresponding stress intensity factors (Tables 25,26,27,28 and 
Figures 11-13). The strain energy release rate, however is as 
consistent as the shear stress, since shear is the dominant stress. 
Therefore G increases as Hi increases m all three cases, hi > h^, 
^ = hj| and hi < h^ . 
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The evaluation of the infinite integrals and consequently the 
numerical solutions of the integral equations as mentioned before, 
becomes very difficult as the thicknesses decrease, resulting in 
costlier computations. For these thin ge cme tries convergence becomes 
very slow, hence no results are given. However, it should be 
mentioned that in this case the stresses tend to oscillate near the 
crack tip, which happens to be in agreement with [ 16 ]. 
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TABLES 
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Table 1. The effect of adhesive thickness on the 
envelope function R1(«) (equation 66a). 

(h.j/a)=(hn/a)=1.0 , a=1 .0 in. 


oc 

(h2/a)=0.0025 

R1(«) 

(hg/aJrO.OOS 

(h2/a)=0.05 

5.878 

2.272879 

2.213900 

1.279204 

10.755 

2.223983 

2.116758 

.763813 

20.510 

2.126674 

1.926407 

.430288 

30.265 

2.030539 

1.744085 

.34 82 93 

40 . 020 

1.935993 

1.572448 

.2 6 5 924 

54.653 

1.797995 

1.339239 

.138605 

64.408 

1.709048 

, 1.201911 

.078680 

74.163 

1.622896 

1.07 8744 

.041523 

83.91 8 

1.539796 

.970431 

.020802 

93.673 

1.459961 

. 876033 

.010035 

103.429 

1.383560 

.7 946 02 

.004705 

113.184 

1.310719 

.724 990 

.00215 8 

122.939 

1.241522 

.6 6 5 952 

.000973 

132.694 

1. 176016 

.616218 

.000432 

142.449 

1.114207 

.574549 

.000189 

152.204 

1.056 074 

.539778 

.000082 

161.959 

1.001 560 

.5 1 0833 

.000035 

171.714 

.950589 

.4 86748 

.000015 

181.46 9 

.903059 

.466669 

.000006 

191.224 

.858851 

.44984 9 

.000003 

200.980 

.817 834 

.435643 

.000001 

210.735 

.779864 

.423500 

.0 

220.4 90 

.744790 

.412957 

.0 

230.245 

.712455 

.403622 

.0 

240.000 

.6 827 01 

.395173 

.0 
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Table 2. The effect of adhesive thickness on the 
envelope function T2(CC) (equation 66d). 

(h^/a)s(hjj/a)=1.0 , a=1.0 in. 


a 

(hg/a) =0.0025 

T2(CC) 

(hg/aJrO.OOS 

(h2/a)=0.05 

5.878 

-2.32256 8 

-2.313114 

-2.133485 

10.755 

-2.314808 

-2.2 97474 

-1.9381 80 

20.510 

-2.2 990 89 

-2.265769 

-1.486 075 

30.265 

-2.283301 

-2.233475 

-1.035337 

40.020 

-2.267406 

-2.200334 

-.654195 

54.653 

-2.2 432 82 

-2.14 8607 

-.2772 56 

64.408 

-2.226963 

-2.112594 

-.143009 

74.163 

-2.210420 

-2 . 0752 91 

-.070119 

83.91 8 

-2.193 6 2 8 

-2.036698 

-.033151 

93 . 67 3 

-2.176563 

-1.996856 

-.015263 

103.429 

-2.159209 

-1.95 5 842 

-.006888 

113-184 

-2.141551 

-1. 913753 

-.003060 

122.93 9 

-2.12357 8 

-1.87 0701 

-.001342 

132.694 

-2.105284 

-1.826807 

-.000582 

142.449 

-2.086666 

-1.782 1 91 

-.000250 

152.204 

-2.067723 

-1.736974 

-.000107 

161.959 

-2.048458 

-1.691273 

-.000045 

171.714 

-2.028877 

-1.645198 

-.00001 9 

181 .46 9 

-2.008986 

-1.5 9 8 856 

-.000008 

191.224 

-1.988795 

-1.552346 

-.000003 

200.980 

-1.96 8316 

-1.505762 

-.000001 

210.735 

-1.947559 

-1.45 91 92 

-.000001 

220.4 90 

-1.926539 

-1.412722 

.0 

230.245 

-1 . 905270 

-1.366429 

.0 

240.000 

-1.883765 

-1.320391 

.0 
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Table 3. The effect of (hg/a) ratio on the adhesive stresses 
<r and t with equal thickness adherends. 

(h 1 /a)=(h 1| /a)=0.35 , a=1.0 in. 


h 2 /a — 
x/a 

0.0025 

-*/<r 0 

0.004 

-*/*o 

0.005 

0.006 

-*/* o 

.99997 

1.799 

1.782 

1.775 

1.767 

.99883 

.363 

.335 

.324 

.316 

.99533 

.263 

.222 

.209 

.200 

.98951 

.246 

.194 

.174 

.162 

.981 37 

.223 

.188 

.167 

.151 

. 97094 

.203 

.176 

.161 

.146 

.95825 

.188 

.162 

.151 

.140 

.94331 

.175 

.148 

.140 

.132 

.92617 

.164 

.135 

.128 

.122 

. 906 87 

. .151 

.123 

.116 

.111 

.88546 

.138 

.112 

.105 

.101 

.86197 

.124 

.102 

.0952 

.0916 

. 83647 

.108 

.0922 

.0859 

.0826 

.80 902 

.0928 

.0827 

.0772 

.0741 

.77967 

.0775 

.0734 

.0690 

.0663 

.74851 

.0631 

.0644 

.0613 

.05 90 

.71560 

.0500 

.0558 

.0539 

.0522 

. 6 81 02 

.0385 

.0 476 

.0470 

.0459 

.64484 

.0290 

.0400 

.0406 

.0401 

.60716 

.0212 

.0331 

.0347 

.034 8 

.56 806 

.0152 

.0270 

.0293 

.0299 

.52764 


.0217 

.0245 

.0255 

.4 85 98 

.00716 

.0171 

.0203 

.0216 

.44319 

.00477 

.0133 

.0166 

.01 81 

' .35460 


.00772 

.0107 

.0123 

.30 902 

.00118 

.00 574 

.00842 

.00995 

.26271 

.00073 

.00421 

.00652 

.007 90 

.16836 

.00021 

.00207 

.00356 

.00453 

.07243 

.0 

.00076 

.00139 

.00182 

.02416 

.0 

.00025 

.00046 

.00060 
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Table 3 . ( cont . ) 


hg/a — ' 

x/a 

0.007 

-*'*o 

0.008 

-*'*o 

0.009 

-*/*o 

for all 
hg/a ratios 

. 99997 

1.759 

1.751 

1.743 

0.0 

. 99883 

.310 

.305 

.301 

0.0 

.99533 

.1 92 

.186 

.180 

0.0 

.98951 

.154 

.148 

.143 

0.0 

.98137 

.140 

.132 

.125 

0.0 

. 970 94 

.135 

.126 

.118 

0.0 

.95825 

.130 

.121 

.114 

0.0 

.94331 

.123 

.116 

.109 

0.0 

.92617 

.116 

.110 

.104 

0.0 

.906 87 

.107 

.102 

.0977 

0.0 

.88546 

.0982 

.0948 

.0912 

0.0 

.86197 

.0892 

.0868 

.0842 

0.0 

.83647 

.0806 

.0789 

.0771 

0.0 

.80 902 

.0724 

.0712 

.0700 

0.0 

.77967 

.0648 

.0638 

.0630 

0.0 

.74851 

.0 577 

.0 570 

.0565 

0.0 

.71560 

.0512 

.0506 

.0503 

0.0 

.6 8102 

.0451 

.0447 

.0446 

0.0 

.64484 

• 0396 

.0393 

.0393 

0.0 

.60716 

.0346 

.0344 

.0345 

0.0 

.5 6 806 

.0300 

.0300 

.0301 

0.0 

.52764 

.0258 

.0259 

.0261 

0.0 

.48598 

.0221 

.0223 

.0225 

0.0 

.44319 

.0187 

.0190 

.0193 

0.0 

.35460 

.0131 

.0135 

.0138 

0.0 

.30 902 

.0107 

.0111 

.0114 

0.0 

.26271 

.00863 

.00902 

.00929 

0.0 

.16836 

.00506 

.00536 

.00555 

0.0 

.07243 

.00207 

.00221 

.00229 

0.0 

.02416 

.00068 

.00073 

.00076 

0.0 
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Table 4 


. The effect of (hg/a) ratio on the adhesive stresses 
and T with equal thickness adherends. 

(hj/a)=(hj|/a)=0.25 , a=1.0 in. 


hg/a — •» 
x/a 

0.0025 

0.004 

-*/*o 

0.005 

0.006 

-T/<T 0 

.99997 

1.515 

1.496 

1.488 

1.480 

.99883 

.305 

.281 

.272 

.265 

.99533 

.221 

.185 

.175 

.167 

.98951 

.205 

.162 

.145 

.135 

.981 37 

.185 

.156 

.138 

.125 

. 97094 

.167 

.145 

.132 

.120 

• 95825 

.153 

.132 

.123 

.114 

.94331 

.141 

.119 

.112 

.106 

.92617 

.130 

.107 

.101 

.0969 

, . 906 87 

.118 

.0960 

.0 908 

.0875 

’ .88546 

.105 

.0860 

.0809 

.0782 

.86197 

.0909 

.0765 

.0718 

.0694 

. 83647 

.0768 

.0675 

.0634 

.0613 

.80 902 

.0629 

.0589 

.0556 

.0538 

.77967 

.0497 

.0506 

.0483 

.0469 

.74851 

.0379 

.0428 

.0416 

.0406 

.71560 

.0278 

.0356 

.0355 

.0350 

.68102 

.01 96 

.02 91 

.0299 

.0298 

.64484 

.0134 

.0233 

.0249 

.0252 

.60716 

.0 0 872 

.0183 

.0204 

.0212 

.56 806 

.00550 

.0142 

.0166 

.0176 

.52764 

.00334 

.0107 

.0133 

.0145 

.485 98 

.001 93 

.0796 

.0105 

.0118 

.44319 

.00110 

.005 81 

.00822 

.00957 

. .35460 

.00029 

.00293 

.00482 

.00605 

.30 902 

.00014 

.00203 

.00362 

. 00 473 

.26271 

.0 

.00138 

.0026 8 

.00363 

.16836 

.0 

.00059 

.00135 

.00197 

.07243 

.0 

.00020 

.00050 

.00077 

.02416 

.0 

.0 

.00016 

.00025 


60 



Table 4 . ( cont . ) 


for all 

i ^/ a — 0.007 0.008 0.009 hg/a ratios 


x/a 


-*'*o 

-*/*o 

v/v 

. 99997 

1.471 

1.463 

1.455 

0.0 

. 99883 

.259 

.255 

.251 

0.0 

.99533 

.160 

.155 

.150 

0.0 

.98951 

.128 

.123 

.119 

0.0 

. 981 37 

.116 

.109 

.103 

0.0 

. 970 94 

.111 

.103 

.0970 

0.0 

.9 5 825 

.106 

.0985 

.0 926 

0.0 

.94331 

.0995 

.0933 

.0879 

0.0 

. 926 17 

.0922 

.0874 

.0827 

0.0 

.906 87 

.0842 

.0807 

.0770 

0.0 

.88546 

.0760 

.0735 

.0709 

0.0 

.86197 

.0678 

.0662 

.0644 

0.0 

.83647 

.0601 

.0591 

.057 9 

0.0 

.80 902 

.0528 

.0522 

.0515 

0.0 

.77967 

.0462 

.0458 

.0455 

0.0 

.74851 

.0401 

.0399 

.0398 

0.0 

.71560 

.0346 

.0346 

.0346 

0.0 

.68102 

.02 97 

.0298 

.0300 

0.0 

.64484 

.0254 

.0255 

.0258 

0.0 

.60716 

.0215 

.0217 

.0220 

0.0 

.56 806 

.0181 

.0184 

.0187 

0.0 

.52764 

.0151 

.0155 

.0158 

0.0 

.48598 

.0125 

.0130 

.0133 

0.0 

.44319 

.0103 

.0108 

.0111 

0.0 

.35460 

.00679 

.00724 

.00758 

0.0 

.30 902 

.00541 

.00584 

.00615 

0.0 

.26271 

.00424 

.00462 

.00490 

0.0 

.16836 

.00238 

.00265 

.00283 

0.0 

.07243 

.00095 

.00107 

.00115 

0.0 

.02416 

.00031 

.00035 

.00038 

0.0 
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Table 5. The effect of (hj/a , hjj/a) ratios on the adhesive 
stresses and 'l with equal thickness adherends. 


x/a 


.99997 
. 99883 
.99533 
.98951 
.981 37 
. 970 94 
.95 825 
. 94331 
.92617 
.90687 
.88546 
.86197 
.83647 
.80 902 
.77967 
.74851 
.71560 
.68102 
.64484 
.60716 
.56 806 
.52764 
.4 85 98 
.44319 
.35460 
.30 902 
.26271 
.16 836 
.07243 
.02416 


(hg/a) =0.0025 , h^hjj , 


0.25 

-*/*o 

0.30 

-*/*o 

0.35 

-T/<T 0 

1.515 

1.663 

1.799 

.305 

.336 

.363 

.221 

.243 

.263 

.205 

.226 

.246 

.185 

.205 

.223 

.167 

.186 

.203 

.153 

.171 

.188 

.141 

.159 

.175 

.130 

.147 

.164 

.118 

.135 

.151 

.105 

.122 

.138 

.0909 

.108 

.124 

.0768 

.0 932 

.108 

.0629 

.0783 

.0928 

.0497 

.0639 

.0775 

.0379 

.0506 

.0631 

.0278 

.0388 

.0500 

.0196 

.0289 

.0385 

.0134 

.0208 

.0290 

.00 872 

.0145 

.0212 

.00550 

.00982 

.0152 

.00334 

.00646 

.0106 

.00193 

.0040 8 

.00716 

.00110 

.00254 

.00477 

.00029 

.00086 

.00192 

.ootm 

.00048 

.00118 

.0 

.00029 

.00073 

.0 

.0 

.00021 

.0 

.0 

.0 

.0 

.0 

.0 


=1.0 in. 

for all 

0.50 1.0 h.j/a ratios 

*/* 0 --r/<r 0 <r/<r 0 


157 

3.059 

0.0 

436 

.619 

0.0 

317 

.451 

0.0 

297 

.424 

0.0 

271 

.389 

0.0 

248 

.358 

0.0 

231 

.337 

0.0 

217 

.321 

0.0 

205 

.307 

0.0 

193 

.2 93 

0.0 

179 

.27 8 

0.0 

164 

.261 

0.0 

148 

.243 

0.0 

131 

.223 

0.0 

114 

.203 

0.0 

0S73 

.182 

0.0 

0817 

.163 

0.0 

0674 

.144 

0.0 

0547 

.127 

0.0 

0437 

.111 

0.0 

0345 

.0964 

0.0 

0268 

.0836 

0.0 

0205 

.0720 

0.0 

0155 

.0618 

0.0 

0 0 845 

.0442 

0.0 

00611 

.0367 

0.0 

00437 

.0300 

0.0 

00200 

.0179 

0.0 

00074 

.00749 

0.0 

00025 

.00250 

0.0 
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Table 6. The effect of (h 2 /a) ratio on the normalized 
stress intensity factors (k-j/k 0 , k 2 /k 0 ) and 
the strain energy release rate, G/G 0 , with 
equal thickness adherends. 

OyaMiyaJsO.BS , k 0 ^ 0 YT 

G^ff^a/Eg , a= 1 . 0 in. 


h 2 /a 0.0025 0.004 0.005 0.006 



k ^/ k 0 0.0 0.0 0.0 

k 2 /k 0 -.01363 -.01356 -.01350 

G/G q .490E-3 .485E-3 .481E-3 
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Table 7 


. The effect of (hg/a) ratio on the normalized 
stress intensity factors (k.j/k 0 , k 2 /k Q ) and 
the strain energy release rate, G/G QJ with 
equal thickness adherends. 

(h^/a)=(h|j/a)=0.25 , k 0 =ff^/T 



G o=*§ a/E 2 

, a= 1 . 0 in. 



h 2 /a 

0.0025 

0.004 

0.005 

0.006 

k l/k 0 

0.0 

0.0 

0.0 

0.0 

k 2 /k 0 

-.01173 

-.01159 

-.01153 

-.01146 

G/Gq 

.363E-3 

•354B-3 

.351 E-3 

.347 E-3 

hg/a 

0.007 

0.008 

0.009 


k 1 /k o 

0.0 

0.0 

0.0 


k 2 /k o 

-.01140 

-.01133 

-.01127 


G/G c 

•343E-3 

•339E-3 

•335E-3 
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Table 8. The effect of (h^/a , h^/a) ratios on the normalized 
stress intensity factors (k.|/k 0 , k 2 /k 0 ) and the. 
strain energy release rate, G/G 0 , with equal 
thickness adherends. 

(h2/a)=0.0025 , k 0 =V^/a , h 1 = h 4 

G 0 =<r^a/E2 , a=1.0 in. 


h.,/a 

0.25 

0.30 

0.35 

0.50 

1.0 

k 1 /k o 

0.0 

0.0 

0.0 

0.0 

0.0 

V k o 

-.01173 

-.01288 

-.01394 

-.01671 

-.02369 

G/G 0 

.363E-3 

.438E-3 

. 51 3E-3 

.737E-3 

. 1 48E-2 
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Table 9. The effect of H , the lumber of unknown functions 
selected in the solution of the system of integral 
equations (Section 3.1) , on the stress intensity 
factor and the strain energy release rate for thin 
geometries and the comparision of (G/G Q ) with the 
plate solution in Ref. [4] . 
hi = h|j=0. 1 25 in., h2=0.0025 in. 

k o =<r oV^~» ^ 0^2 ‘ 

G plate=°- 2581 ’ * 0 = 2E+it lb/in2 * 




N=26 

N=29 

N=33 

- Plate 
Soln. 

asl.O 

k 1 /k o 

k 2 /k o 

0.0 

-.00816 

0.0 

-.0 0 821 

0.0 

-.00818 

O 

• 

o 


G/Gq 

. 176 E-3 

. 178E-3 

. 177E-3 

.181 E-3 



N=37 

N=40 

N=45 



k 1 /k o 

V k o 

0.0 

-.0 0 57 8 

0.0 

-.00581 

0.0 

-.0 0 57 9 

0.0 

) 

■ 

G/Gq 

. 8 82 E-4 

. 8 91 E-4 

. 8 85 E-4 

| . 90 3 E-4 

1 

1 

l 



N=4l 

N=45 



a=3.0 

k 1 /k o 

k 2 /k o 

0.0 

-.00465 

0.0 

-.00471 


0.0 


g/g 0 

.57 IE-4 

.5 85 E-4 


.602E-4 
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Table 10. The effect of bond length (a), on the adhesive 
stresses <T and T , with equal thickness 
adherends, for thin geometries. 

h.|=hjj=0. 125 in. , ^=0.0025 in. 


a — - 
x/a 

1.0 

-*/*o 

2.0 

-T/ff 0 

3.0 

1.0, 2. 0,3.0 

.99997 

1.052 

.748 

.609 

0.0 

.99938 

.265 

.208 

.189 

0.0 

.99751 

.172 

.151 

.144 

0.0 

. 99440 

.144 

.137 

.124 

0.0 

.99005 

.138 

.121 

.110 

0.0 

.98447 

.131 

.106 

.0990 

0.0 

.97766 

.121 

.0943 

.0861 

0.0 

. 96 963 

.109 

.0832 

.0704 

0.0 

.96039 

.097 9 

.0718 

.0526 

0.0 

.94996 

.0872 

.06 00 

.0350 

0.0 

.93834 

.0774 

.0479 

.0202 

0.0 

. 92556 

.0684 

.0362 

.00993 

0.0 

.91 162 

.0600 

.0258 

.00408 

0.0 

.896 55 

.0522 

.0170 

.00139 

0.0 

.88036 

.0448 

.0104 

.00038 

0.0 

.86307 

.0380 

.00 578 

.0 

0.0 

. 8447 1 

.0317 

.00293 

.0 

0.0 

.82529 

.0261 

.00135 

.0 

0.0 

.80485 

.0210 

.00054 

.0 

0.0 

.7 8340 

.0167 

.00020 

.0 

0.0 

.76098 

.0129 

.0 

.0 

0.0 

.73761 

.00983 

.0 

.0 

0.0 

' .71332 

.007 31 

.0 

.0 

0.0 

.66211 

.00378 

.0 

.0 

0.0 

.60759 

.00178 

.0 

.0 

0.0 

.55006 

.00076 

.0 

.0 

0.0 

.4 8978 

.00029 

.0 

.0 

0.0 

.42706 

.00010 

.0 

.0 

0.0 

.39488 

.0 

.0 

.0 

0.0 

.01765 

.0 

.0 

.0 

0.0 
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Table 11. The effect of (hg/a) ratio on the adhesive 
stresses 9 and T with different thickness 
adherends. 


(h-,/a)=1.0 , (hjj/a)=0.5 , a=1.0 in. 


(h2/a)=0.03 


-t/9 . 


9 / 9 . 


(hg/aJsO.Oil 


- 1 / 9 , 


9 / 9 , 


x/a 


.99997 
.99883 
.99533 
.98951 
.98137 
. 970 94 
.95825 
.94331 
.92617 
. 906 87 
.88546 
.86197 
.83647 
.80 902 
.77967 
.74851 
.71560 
.68102 
.64484 
.60716 
.56806 
.52764 
.4 85 98 
.44319 
.35460 
.30 902 
.26271 
.16836 
.07243 
.02416 


2.174 
.356 
.190 
.139 
.116 
.101 
.0912 
.0844 
.0801 
.0771 
.0745 
.0717 
.0689 
.065 9 
.0629 
.05 97 
.0565 
.0532 
.0499 
.0466 
.0432 
.0398 
.0364 
.0330 
.0261 
.0226 
.0191 
.0122 
.00523 
.00174 


-.523 

-.0891 

-.0527 

-.0444 

-.0426 

-.0420 

-.0409 

-.0389 

-.0363 

-.0334 

-.0301 

-.0267 

-.0232 

-0197 

-.0162 

-.0128 

-.00950 

-.00632 

-.00326 

-.00034 

.00243 

.00505 

.00750 

.00977 

.0137 

.0154 

.0168 

.0190 

.0203 

.0205 


2.023 
.330 
.173 
.124 
.101 
.0879 
.0786 
.0716 
.0666 
.0631 
.0606 
.05 83 
.0561 
.0538 
.0514 
.0489 
.0463 
.0437 
.041 1 
.0384 
.0356 
.0329 
.0301 
.0273 
.0216 
.0188 
.0159 
.0102 
.00436 
.00145 


-.464 
-.077 8 
-.0441 
-.0354 
-.0327 
-.0318 
-.0312 
-.0301 
-.0284 
-.0264 
-.0240 
-.0215 
-.0189 
-.0162 
-.0134 
-.0107 
-.007 94 
-.00531 
-.00276 
-.00033 
.00197 
.00412 
.00612 
.00795 
.0111 
.0124 
.0135 
.0152 
.0162 
.0163 
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Table 11 . ( eont . ) 


x/a 

( h 2 / a ) 

= 0.05 

*'*o 

( h 2 / a ) 

-*/<r o 

= 0.06 

<r/<r 0 

.99997 

1.905 

-.41 9 

1.806 

-.385 

. 99883 

.309 

-.0695 

.292 

-.0635 

.99533 

.160 

-.0384 

.151 

-.0344 

.98951 

.113 

-.02 97 

.105 

-.0260 

.98137 

.0912 

-.0266 

.0838 

-.0227 

. 970 94 

.0786 

-.0255 

.0716 

-.0214 

.95825 

.0700 

-.0249 

.0635 

-.0207 

.94331 

.0635 

-.0242 

.0574 

-.0201 

. 926 17 

.0583 

-.0231 

.0525 

-.0193 

.90687 

.0545 

-.0216 

.0485 

-.0182 

.88546 

.0516 

-.0199 

.0455 

-.0169 

.86197 

.0494 

-.0180 

.0431 

-.01 53 

.83647 

.0475 

-.0160 

.0412 

-.0137 

.80 902 

.0455 

-.0138 

.0394 

-.0119 

.77967 

.0436 

-.0116 

.0377 

-.0101 

.74851 

.0415 

-.00940 

.0360 

-.0 0 822 

.71560 

.0394 

-.00716 

.0342 

-.00633 

.68102 

.0372 

-.004 96 

.0323 

-.00445 

.64484 

.0350 

-.00280 

.0304 

-.00260 

.60716 

.0327 

-.00072 

.0285 

-.00080 

.56 806 

.0304 

.00126 

.0265 

.00092 

. 52764 

.02 81 

.00314 

.0245 

. 002 56 

.4 84 9 8 

.0257 

.00490 

.0224 

.00410 

.44319 

.0233 

.00652 

.0204 

.00553 

.35460 

.0185 

.00933 

.0162 

.0 0 800 

.30 902 

.0160 

.0105 

.0141 

.00 903 

.26271 

.0136 

.0115 

.0119 

.00993 

.16 836 

.00866 

.0130 

.00762 

.0113 

.07243 

.00372 

.0139 

.00327 

.0121 

.02416 

.00124 

.0141 

.00109 

.0122 
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Table 11 . ( cont . ) 


( h 2 / a )= 0.07 ( h 2 / a )= 0.1 

-*/*o <r/<T 0 -T/<r c «7* 0 


x/a 


.99997 
. 99883 
.99533 
.98951 
.98137 
.97094 
.95825 
.94331 
.92617 
.906 87 
.88546 
.86197 
.83647 
.80 902 
.77967 
.74851 
.71560 
.68102 
.64484 
.60716 
.56 806 
.52764 
.48598 
.44319 
.35460 
.30902 
.26271 
.16 836 
.07243 
.02416 


1.723 
.278 
.143 
.0990 
.0781 
.0662 
.0585 
.0527 
.0481 
.0443 
.0411 
.0386 
.0366 
.034 9 
.0333 
.0318 
.0302 
.0286 
.026 9 
.0252 
.0235 
.0217 
.0199 
.0180 
.0143 
.0124 
.0105 
.00673 
.00289 
.000 96 


-.357 
-.05 86 
-.0314 
-.0233 
-.0199 
-.0184 
-.0176 
-.0171 
-.0164 
-.0156 
-.0146 
-.0133 
-.0120 
-.0105 
-.00895 
-.00736 
-.0 0 57 5 
-.00412 
-.00251 
-.000 93 
.00060 
.00206 
.00345 
.00474 
.00700 
.00796 
.00879 
.0101 
.0108 
. 0109 ' 


1.535 
.247 
.126 
.085 9 
.0666 
.0555 
.0483 
.0432 
.0392 
.035 9 
.0330 
.0305 
.0284 
.0266 
.0250 
.0237 
.0224 
.0211 
.0199 
.0187 
.0174 
.0161 
.0147 
.0134 
.0107 
.00926 
.00785 
.00501 
.00215 
.00072 


-.299 
-.04 86 
-.0254 
-.0182 
-.0149 
-.0132 
-.0123 
-.0117 
-.0112 
-.0106 
-.0100 
-.00 927 
-.00842 
-.00746 
-.00644 
-.00537 
-.00426 
-.00313 
-.00199 
-.00087 
.00024 
.00131 
.00233 
.00330 
.00501 
.00 575 
.00639 
.00737 
.007 93 
.00804 
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Table 12. The effect of bond length (a), on the adhesive 
stresses CP and f , with different thickness 
adherends. 


hjsl.O in. , hj|=0.5 in. , hgsO.I in. 


a=2.0 


-^/<T 0 */* Q 


a=3. 0 

<r/cr. 


a=4.0 

-•T/C„ <F /<£, 


x/a 


.99997 
.99883 
.99533 
.98951 
.981 37 
.970 94 
.95825 
.94331 
.92617 
. 90687 
.88546 
.86197 
.83647 
'. 80 902 
.77967 
.74851 
.71560 
.6 8102 
.64484 
.60716 
.56 806 
. 52764 
.4 85 9 8 
.44319 
.35460 
.30 902 
.26271 
.16836 
.07243 
.02416 


1.511 
.245 
.127 
.0894 
.0718 
.0616 
.0547 
.04 93 
.0451 
.0419 
.0395 
.0376 
.0359 
.0343 
.0326 
.0310 
.0292 
.0275 
.0257 
.0240 
.0222 
.0204 
.0186 
.0168 
.0133 
.0115 
.00971 
.00617 
.00264 
.00088 


-.414 
-.06 86 
-.0377 
-.0289 
-.0256 
-.0242 
-.0233 
-.0222 
-.0207 
-.0189 
-.0168 
-.0145 
-.0121 
-.00 961 
-.00720 
-.004 86 
-.00265 
-.00061 
.00122 
.002 82 
.0041 8 
.00530 
.00619 
.00686 
.00766 
.00784 
.00793 
.00792 
.00784 
.00782 


1.333 
.218 
.115 
.0827 
.0679 
.0588 
.0522 
.0474 
.0441 
.0417 
.0397 
.0378 
.0358 
.0338 
.0317 
.02 97 
.0277 
.02 57 
.0237 
.0217 
.0198 
.0179 
.0161 
.0143 
.0110 
.00 935 
.00780 
.004 86 
.00205 
.00068 


-.375 
-.0632 
-.0362 
-.02 94 
-.0271 
-.0259 
-.0243 
-.0221 
-.0194 
-.0163 
-.0131 
-.00 979 
-.00660 
-.00361 
-.00092 
.00139 
.00329 
.00476 
.0 0 57 9 
.00643 
.00670 
.00669 
.00643 
.00602 
.00492 
.00435 
.00381 
.00295 
.00244 
.00234 


1.176 
.193 
.103 
.0761 
.0631 
.0547 
.0487 
.044 9 
.0423 
.0400 
.0377 
.0354 
.0331 
.0308 
.0285 
.0262 
.0240 
.0218 
.0197 
.0177 
.0157 
.0139 
.0122 
.0106 
.00771 
.00645 
.00528 
.00319 
.00133 
.000 439 


-.333 
-.0570 
-.0340 
-.0288 
-.0271 
-.0254 
-.0229 
-.0197 
-.0160 
-.0122 
-.0 0 834 
-.00 472 
-.00147 
.00129 
.00348 
.00507 
.00609 
.0065 8 
.00661 
.0062 8 
.0 0570 
.00498 
.00420 
.00344 
.00220 
.00177 
.00145 
.00112 
.00101 
.00099 
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Table 


x/a 


.99997 
. 99883 
.99533 
.98951 
.981 37 
.970 94 
.95 825 
. 94331 
. 926 17 
.906 87 
.88546 
.86197 
>83647 
.80 902 
.77967 
.74851 
.71560 
.68102 
.64484 
.60716 
.56 806 
.52764 
.4 85 98 
.44319 
.35460 
.30 902 
.26271 
.21578 
.16 836 
.072 43 
.02416 


3. The effect of (h^/a) ratio on the adhesive 
stresses (T and t with different thickness 
adherends. 


(h]j/a)=0.5 » 

(h 2 /a)=0. 

1 , a=1 .0 

in. 


(h 1 /a)sl.o 

(h 1 /a)=2.0 

(h^/a)=3.0 

-*/<T 0 

<T/<T 0 

-f /* 0 


-T/<F 0 

<r/<r 0 

1.535 

-.299 

1.605 

-.391 

1.620 

-.405 

.247 

-.04 86 

.259 

-.0637 

.261 

-.065 9 

.126 

-.0254 

.131 

-.0333 

.133 

-.0344 

.085 9 

-.0182 

.0899 

-.0238 

.0 907 

-.02 46 

.0666 

-.0149 

.0698 

-.0195 

.0704 

-.0202 

.0555 

-.0132 

.0582 

-.0173 

.05 87 

-.0179 

.0483 

-.0123 

.0506 

-.0160 

.0511 

-.0166 

.0432 

-.0117 

.0453 

-.0152 

.0457 

-.0158 

.0392 

-.0112 

.0412 

-.0146 

.0415 

-.0151 

.0359 

-.01 06 

.0377 

-.0139 

.0380 

-.0144 

.0330 

-.0100 

.0347 

-.0131 

.0350 

-.0135 

.0305 

-.00927 

.0322 

-.0121 

.0324 

-.0125 

.0284 

-.00842 

.0299 

-.0110 

.0302 

-.0114 

.0266 

-.00746 

.0281 

-.00 971 

.02 83 

-.0101 

.0250 

-.00644 

.0264 

-.00837 

.0267 

-.00868 

.0237 

-.00537 

.0250 

-.006 96 

.02 53 

-.00723 

.0224 

-.00426 

.0237 

-.00551 

.0239 

-.0 0 57 3 

.0211 

-.00313 

.0224 

-.00403 

.0226 

-.00420 

.0199 

-.00199 

.0211 

-.00255 

.0213 

-.00267 

.0187 

-.00087 

.0198 

-.00108 

.0200 

-.00115 

.0174 

.00024 

.0185 

.00035 

.0187 

.00034 

.0161 

.00131 

.0171 

.00174 

.0173 

.00178 

.0147 

.00233 

.0157 

.00307 

.0159 

.00316 

.0134 

.00330 

.0143 

.00432 

.0145 

.00446 

.0107 

.00501 

.0113 

.00653 

.0115 

.00677 

.00 926 

.00 575 

.00980 

.00748 

.0100 

.00776 . 

.00785 

•00639 

.00830 

.00831 

.00850 

.00862 

.00644 

.006 93 

.006 81 

.00 901 

.00698 

.00 936 

.00501 

.00737 

.0 0531 

.00958 

.00544 

.00995 

.00215 

.007 93 

.00229 

.0103 

.00234 

.0107 

.00071 

.00804 

.00076 

.0104 

.0007 8 

.0109 
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Table 


x/a 


.99997 
.99883 
.99533 
.98951 
.98137 
. 97094 
.95825 
.94331 
. 926 17 
.90687 
.88546 
.86197 
.83647 
. 80 902 
.77967 
.74 851 
.71560 
.68102 
.64484 
.60716 
.5 6 806 
. 52764 
.485 98 
.44319 
.35460 
.30 902 
.26271 
.16 836 
.07243 
.02416 


14. The effect of (hg/a) ratio on the adhesive 
stresses 9 and * with different thickness 
adherents. 


(h.,/a)=0.5 , 

(h 4 /a)=1. 

0 , a=1.0 

in. 


uya) 

=0.03 

(hg/a) 

=0.04 

(hg/a) 

=0.05 

-*/*o 

9/9 Q 

-*/*o 

9/9 Q 

-*'*0 

9/9 Q 

2.173 

.525 

2.023 

.464 

1.904 

.419 

.356 

.0895 

.330 

.0778 

• 309 

.06 96 

.190 

.0529 

.173 

.0441 

.160 

.0385 

.139 

.0446 

.124 

.0354 

.113 

.0298 

.115 

.0427 

.101 

.0323 

.0911 

.0267 

.101 

.0421 

.0879 

.0318 

.0785 

.0255 

.0911 

.0410 

.0786 

.0312 

.0700 

.0249 

.0843 

.0390 

.0716 

.0301 

.0634 

.0242 

.0801 

.0364 

.0666 

.0284 

.0583 

.0231 

.0771 

.0334 

.0631 

.0264 

.0544 

.0216 

.0744 

.0301 

.0606 

.0240 

.0516 

.0199 

.0717 

.0266 

.05 83 

.0215 

.04 94 

.0180 

.0688 

.0231 

.0561 

.0189 

.0474 

.0159 

.0659 

.0195 

.0538 

.0162 

.0455 

.0138 

.0628 

.0160 

.0514 

.0134 

.0435 

.0116 

.05 96 

.01 26 

.0489 

.0107 

.0415 

.00 931 

.0 564 

.00931 

.0463 

.00794 

.0393 

.00707 

.0532 

.00612 

.0437 

.00531 

.0372 

.004 86 

.0499 

.00306 

.0411 

.00276 

.0350 

.00271 

.0465 

.00014 

.0384 

.00033 

.0327 

.00064 

.0432 

-.0 0262 

.0356 

-.00197 

.0304 

-.00134 

.0398 

-.00521 

.032 9 

-.00412 

.02 81 

-.00321 

.0364 

-.00763 

.0301 

-.00612 

.0257 

-.00495 

.0330 

-.00985 

.0273 

-.00795 

.0233 

-.00655 

.0261 

-.0137 

.0216 

-.0111 

.0185 

-.00931 

.0227 

-.0153 

.0188 

-.0124 

.0161 

-.0105 

.0192 

-.0167 

.0159 

-.0135 

.0137 

-.0114 

.0123 

-.0188 

.0102 

-.01 52 

.00872 

-.012 9 

.00526 

-.0201 

.00436 

-.0162 

.00374 

-.0138 

.00175 

-.0203 

.00145 

-.0163 

.00125 

-.0140 
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Table 14. (cont. ) 


(h 2 /a)=0.06 




9 / 9 , 


(h 2 /a)=0.07 


- 9 / 9 , 


9 / 9 , 


(h 2 /a)=0. 1 


— 9 / 9 . 


9 / 9 , 


x/a 


.99997 
.997 95 
.99179 
.98156 
.9672 9 
.94 906 
.92692 
.900 97 
.87132 
.83809 
.80141 
-76145 
.71835 
.67230 
.62349 
.57212 
.51839 
. 46254 
.40478 
.34537 
.28453 
.22252 
.15960 
.0 9602 
.03205 


1.806 

.222 

.117 

.0841 

.0688 

.0595 

.0527 

.0476 

.0440 

.0413 

.0390 

.0367 

.0343 

.0319 

.0293 

.0267 

.0240 

.0213 

.0186 

.0158 

.0130 

.0101 

.00723 

.00434 

.00145 


.385 
.0489 
.0281 
.0228 
.0211 
.0203 
.0193 
.0178 
.0160 
.0138 
.0115 
.00 900 
.00650 
.00401 
.00159 
-.00073 
-.00290 
-.00489 
-.00667 
-.00822 
-.00953 
-.0106 
-.0114 
-.0120 
-.0122 


1.723 

.211 

.110 

.0784 

.0635 

.0547 

.0483 

.0433 

.0395 

.0367 

.0345 

.0324 

.0303 

.0282 

.0259 

.0236 

.0213 

.0189 

.0164 

.0140 

.0115 

.00894 

.00640 

.00384 

.00128 


.358 
.0450 
.0253 
.0200 
.0181 
.0173 
.0165 
.01 53 
.0138 
.0120 
.0100 
.007S7 
.00583 
.00367 
.00155 
-.00049 
-.00242 
-.00419 
-.00579 
-.00719 
-.0 0 83 8 
-.00 935 
-.0101 
-.0106 
-.0108 


1.536 

.187 

.0962 

.0670 

.0530 

.0449 
.0394 
.0351 
.0315 
.0285 
.0262 
.0242 
.0225 
.0209 
.0192 
.0175 
.0158 
.0140 
.0122 
.0103 
.00842 
.00656 
.00470 
.002 83 
.00094 


.2 99 
.0371 
.0200 
.0150 
.0129 
.0119 
.0112 
.0105 
.00956 
.00845 
.00718 
.00 578 
.0 0432 
.00282 
.00132 
-.00016 
-.00157 
-.00289 
-.00409 
-.00516 
-.00608 
-.00683 
-.00741 
-.00780 
-.00799 



Table 15. The effect of (hjj/a) ratio on the adhesive 
stresses G and with different thickness 
adherends. 

(h^/a)=0.5 , (h 2 /a)=0.07 , a=1.0 in. 


(h 2} /a)=1.0 

-*T/G G/G 

o o 


<iya)=2.0 


-■T/tf 


G/G. 


(h 4 /a)=3.0 


- t/G . 


G/G, 


x/a 


. 99997 
.99795 
.99179 
.98156 
• 96729 
.94 906 
.92692 
.900 97 
.87132 
_. 83809 
.80141 
.76145 
.71835 
.67230 
.62349 
.57212 
.51839 
.46254 
.40478 
.34537 
.2 8453 
.22252 
.15960 
.0 9602 
.03205 


1.723 
.211 
.110 
.0784 
.0635 
.0547 
.0483 
.0433 
.0395 
.0367 
.0345 
.0324 
.0303 
.02 82 
.0259 
.0236 
.0213 
.0189 
.0164 
.0140 
.0115 
.00894 
.00640 
.00384 
.00128 


.358 
.0450 
.0253 
.0200 
.0181 
.0173 
.0165 
.01 53 
.0138 
.0120 
.0100 
.00797 
.00583 
.00367 
.00155 
-.00049 
-.00242 
-.00419 
-.00579 
-.00719 
-.00838 
-.00 935 
-.0101 
-.0106 
-.0108 


1.81 8 
.223 
.116 
.0829 
.0672 
.0579 
.0512 
.0460 
.0420 
.0391 
.0368 
.0346 
.0325 
.0302 
.027 9 
.0254 
.0230 
.0204 
.0178 
.0151 
.0124 
.00 972 
.00696 
.00419 
.00140 


.47 9 
.06 03 
.0339 
.0268 
.0243 
.0231 
.0220 
.0205 
.0185 
.0161 
.0134 
.0106 
.007 7 8 
.00489 
.00206 
-.00066 
-.00323 
-.00 560 
-.00774 
-.00 961 
-.0112 
-.0125 
-.0135 
-.0141 
-.0145 


1 . 840 
.226 
.118 
.0839 
.0680 
.05 86 
.0518 
.0465 
.0426 
.0396 
.0372 
.0351 
.0329 
.0306 
.0282 
.0258 
.0233 
.0207 
.0180 
.0154 
.0126 
.00986 
.00707 
.00425 
.00142 


.4 99 
.0628 
.0352 
.0278 
.0252 
.0241 
.0229 
.0213 
.0192 
.0167 
.0140 
.0110 
.00806 
.00506 
.00212 
-.00071 
-.00337 
-.005 83 
-.0 0 804 
-.00999 
-.0116 
-.0130 
-.0140 
-.0147 
-.0150 
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Table 16. The effect of (h 2 /a) ratio on the normalized 
stress intensity factors (k^/k 0 , k 2 /k Q ) and 
the strain energy release rate , G/G of with 
different thickness adherends. 

(h^/aJsl.O , ( h^/a ) =0. 5 , 

G o =e ’o a/,E 2 * 33 in * 


h 2 /a 

0.03 

0.04 

0.05 

0.06 

k l/k 0 

-.00405 

-.00359 

-.00324 

-.00298 

k 2 /k 0 

-.0168 

-.0157 

-.0148 

-.0140 

G/G 0 

.745E-3 

.650E-3 

.578E-3 

•517E-3 


i^/a 

0.07 

0. 1 

k 1 /k o 

-.00277 

-.00232 

k 2 /k 0 

-.0133 

-.0119 

G/Gq 

.467E-3 

.374E-3 
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Table 17. The effect of bond length (a) , on the normalized 
stress intensity factors (k«/k Q , k 2 /k 0 ) and the 
strain energy release rate, G/G 0 , with 
different thickness adherends. 

h^sl.O in., l^sO.I in., hj,=0.5 in., k 0 =ff^/a» 
G 0=*0 a/E 2 • 



a=2.0 

a=3.0 

a=4.0 

V k o 

-.00321 

-.002 91 

-.0025 8 

k 2 / k 0 

-.0117 

i 

• 

o 
— * 
o 

UJ 

-.0091 1 

G/G 0 

.36 IB-3 

• 2 80E-3 

.219E-3 


Table 18. The effect of (h^/a) ratio on the normalized 
stress intensity factors (k^k^ k 2 /k 0 ) and 
the strain energy release rate , G/G Q , with 
different thickness adherends. 

(h^/a)=0.5 , (hg/aJsO.I , k Q =ff 0 a 

G 0 =ff|a/E 2 f a=1 . 0 in. 


h^a 

1.0 

2.0 

3.0 

k 1 /k o 

-.00232 

-.00303 

-.00314 

k 2 /k o 

-.0119 

-.0124 

-.0125 

G/G q 

.374E-3 

.406E-3 

.41 2E-3 
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Table 19. The effect of (h 2 /a) ratio on the normalized 
stress intensity factors (k.j/k 0 , k 2 /k Q ) and 
the strain energy release rate, G/G 0 , with 
different thickness adherends. 



Oya)=0.5 , 

(h 4 /a)=1.0 

* k 0 =<T 0 \/r 



- 

a=1 .0 in. 



h 2 /a 

0.03 

0.04 

0.05 

0.06 


.00407 

.00360 

.00325 

.00298 

k 2 /k 0 

-.0168 

-.0157 

-.0148 

-.0140 

G/G 0 

.789E-3 

.685E-3 

.606E-3 

.541 E-3 


hg/a 

0.07 

0.1 

k i/k 0 

.0 0277 

.00232 

k 2 /k 0 

-.0133 

-.0119 

g/g 0 

.4 87 E-3 

.388E-3 
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Table 20 


. The effect of (hjj/a) ratio on the normalized 
stress intensity factors (k.j/k 0 , k 2 /k Q ) and 
the strain energy release rate, G/G o , with 
different thickness adherends. 

(h.,/a)=0.5 , ( hg/aJsO.O? , k Q =9jfiT 

G 0 =<f|a/E 2 , a=1 .0 in. 


h^/a 

1.0 

2.0 

3.0 

k l/k 0 

.00277 

.00371 

.00386 

k 2 /k 0 

-.0138 

-.0141 

-.0143 

G/G q 

.523E-3 

.56 IE-3 

.579E-3 
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Table 


x/a 


.99997 
.99795 
.99179 
.98156 
.96729 
.94 906 
.92692 
. .900 97 
.87132 
.83 809 
.80141 
.76145 
.71835 
.67230 
.62349 
.57212 
.5183 9 
. 46254 
.40478 
.34537 
.28453 
.22252 
.15960 
.0 9602 
.03205 


21. The effect of the material properties on the 
adhesive stresses <T and with different 
thickness adherends. 

(h 1 /a)s1.0 , (hjj/a)=0.5 , 

(hg/aJsO.OT , a=1.0 in. 


H-,=M 

-*'*o 

4/2 

*1 = 

*4 

<r/<r 0 

P!=2Mi, 

-t/<t 0 <r/<r 0 

1.582 

-.205 

1.723 

-.358 

1.813 

-.459 

.194 

-.0258 

.211 

-.0450 

.222 

-.0579 

.101 

-.0145 

.110 

-.0253 

.116 

-.0326 

.0716 

-.0114 

.0784 

-.0200 

.0828 

-.0258 

.0578 

-.0104 

.0635 

-.0181 

.0672 

-.0234 

.04 96 

-.00990 

.0547 

-.0173 

.05 80 

-.0224 

.0436 

-.0 0 946 

.0483 

-.0165 

.0513 

-.0213 

.0389 

-.00880 

.0433 

-.0153 

.0461 

-.0199 

.0354 

-.00793 

.0395 

-.0138 

.0422 

-.0180 

.0327 

-.00689 

.0367 

-.0120 

.0393 

-.0157 

.0306 

-.00574 

.0345 

-.0100 

.0370 

-.0132 

.02 86 

-.00454 

.0324 

-.00797 

.034 9 

-.0105 

.0267 

-.00331 

.0303 

-.00583 

.0327 

-.00771 

.02 47 

-.00207 

.02 82 

-.00367 

.0305 

-.004 91 

.0226 

-.00086 

.0259 

-.00155 

.0281 

-.00214 

.0206 

.00030 

.0236 

.00049 

.02 57 

.00054 

.0184 

.00140 

.0213 

.00242 

.0232 

.00307 

.0163 

.00241 

.0189 

.00419 

.0206 

.00540 

.0141 

.00332 

.0164 

.00579 

.0180 

.00751 

.0120 

.00411 

.0140 

.00719 

.0153 

.00 937 

.00981 

.0047 9 

.0115 

.0 0 838 

.0126 

.0109 

.00763 

.00534 

.00894 

.00 935 

.00982 

.0122 

.00545 

.00576 

.00640 

.0101 

.007 04 

.0132 

.00327 

.00604 

.00384 

.0106 

.00423 

.0139 

.00109 

.00618 

.00128 

.0108 

.00141 

.0142 
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Table 


x/a 


. 99997 
. 997 95 
.99179 
.98156 
.96729 
.94 906 
. 926 92 
- . 90097 
.87132 
. 83809 
.80141 
.76145 
.71835 
.67230 
.62349 
.57212 
.5183 9 
. 46254 
.40478 
.34537 
.28453 
.22252 
.15960 
.09602 
.03205 


22. The effect of the material properties on the 
adhesive stresses 9 and 7 with different 
thickness adherends. 

{h-|/a)s0.5 , (hjj/a)=1 .0 , 

(hg/aJsO.O? , a=1 .0 in. 


*1= 

-*/* 0 

M 4 /2 

9/9 Q 

*1 = 
-*/<F 0 

*4 

9/9 Q 

*1 = 
-*/<r 0 

2*4 

9/9 q 

1.570 

.534 

1.723 

.358 

1.839 

.172 

.192 

.0670 

.211 

.0450 

.226 

.0216 

.100 

.0374 

.110 

.0253 

.118 

.0122 

.0712 

.02 93 

.0784 

.0200 

.0840 

.00971 

.0575 

.0263 

.0635 

.0181 

.0682 

.00890 

.04 93 

.0247 

.0547 

.0173 

.0588 

.00858 

.0435 

.0233 

.0483 

.0165 

.0520 

.00828 

.0389 

.0213 

.0433 

.01 53 

.0467 

.00780 

.0354 

.0189 

.0395 

.0138 

.0428 

.00712 

.0329 

.0161 

.0367 

.0120 

.0398 

.00628 

.0308 

.0132 

.0345 

.0100 

.0374 

.00534 

.0289 

.0102 

.0324 

.00797 

.0352 

.00433 

.0270 

.007 0 8 

.0303 

.00583 

.0330 

.00326 

.0251 

.00408 

.02 82 

.00367 

.0307 

.00216 

.0231 

.00120 

.0259 

.00155 

.0283 

.00107 

.0210 

-.00150 

.0236 

-.00049 

.025 9 

-.00001 

.0189 

-.00399 

.0213 

-.00242 

.0233 

-.00104 

.0168 

-.00622 

.0189 

-.00419 

.0207 

-.00201 

.0146 

-.0081 8 

.0164 

-.00579 

.0180 

-.00290 

.0124 

-.00986 

.0140 

-.00719 

.01 53 

-.00370 

.0102 

-.0113 

.0115 

-.0 0 83 8 

.0126 

-.00438 

.00797 

-.0124 

.00894 

-.00 935 

.00981 

-.004 94 

.00571 

-.0132 

.00640 

-.0101 

.00702 

-.00537 

.00343 

-.0138 

.00384 

-.0106 

.00422 

-.00 567 

.00115 

-.0141 

.00128 

-.0108 

.00141 

-.00580 
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Table 


x/a 


. 99997 
.99883 
.99533 
.98951 
.98137 
.970 94 
.95825 
. 94331 
.92617 
.90687 
.88546 
.86197 
.83647 
. 80 902 
.77967 
.74 851 
.71560 
.68102 
.64484 
.60716 
.56806 
.52764 
.4 85 98 
.44319 
.35460 
.30 902 
.26271 
. 16836 
.07243 
.02416 


23. The effect of the material properties on the 
adhesive stresses «T and t with equal thickness 
adherends. 

U^/aMhjj/a^O.S , (h 2 /a)=0.07 , a=1.0 in. 


-T/<T 0 <T/tf 0 

*1 s *4 

-*/<T 0 

<S/<S 0 

.368 

.204 

1.529 

0.0 

1.669 

-.198 

.221 

.0335 

.247 

0.0 

.270 

-.0324 

.113 

.0178 

.127 

0.0 

.138 

-.0174 

.0781 

.0131 

.0876 

0.0 

.095 9 

-.0129 

.0614 

.0111 

.0691 

0.0 

.0757 

-.0110 

.0518 

.0102 

.0584 

0.0 

.0642 

-.0102 

.0455 

.00 959 

.0515 

0.0 

.0566 

-.00973 

.0408 

.00 915 

.0463 

0.0 

.0511 

-.00940 

.0370 

.00867 

.0421 

0.0 

.0466 

-.00903 

.0338 

.0 0 808 

.0386 

0.0 

.042 9 

-.00855 

.0312 

.00739 

.0357 

0.0 

.0398 

-.007 95 

.02 91 

.00663 

.0334 

0.0 

.0374 

-.00727 

.0273 

.00581 

.0316 

0.0 

.0354 

-.00650 

.0259 

.004 96 

.0300 

0.0 

.0337 

-.00569 

.0245 

.00409 

.0285 

0.0 

.0321 

-.00484 

.0232 

. 00322 

.0271 

0.0 

.0306 

-.00395 

.0219 

.00235 

.0256 

0.0 

.0290 

-.00306 

.0205 

.00151 

.0241 

0.0 

.0274 

-.00217 

.0192 

.00070 

.0226 

0.0 

.0258 

-.00128 

.0178 

-.00008 

.0211 

0.0 

.0241 

-.00042 

.0164 

-.00080 

.0195 

0.0 

.0224 

.00041 

.0151 

-.00147 

.0179 

0.0 

.0207 

.00121 

.0137 

-.00208 

.0164 

0.0 

.0189 

.00195 

.0124 

-.00264 

.0148 

0.0 

.0171 

.00264 

.00968 

-.00356 

.0116 

0.0 

.0136 

.00384 

.00 836 

-.00393 

.0101 

0.0 

.0118 

.00434 

.00705 

-.00424 

.00851 

0.0 

.00996 

.00477 

.00446 

-.00469 

.00540 

0.0 

.00634 

.00542 

.00191 

-.00494 

.00231 

0.0 

.0 0272 

.0 0 57 8 

.00064 

-.00499 

.00077 

0.0 

.000 91 

.005 8 
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Table 


x/a 


.99997 
.99883 
.9953 3 
.98951 
.98137 
.970 94 
.95 825 
. 94331 
.92617 
..90687 
.88546 
.86197 
.83647 
.80 902 
.77967 
.74851 
.71560 
.6 8102 
.64484 
.60716 
.56806 
. 52764 
.4 85 98 
.44319 
.35460 
.30 902 
.26271 
.16836 
.07243 
.02416 


24. The effect of the material properties on the 
adhesive stresses and * with equal thickness 
adherends. 

(h^aJsOty/aJsl.O , (h 2 /a)=0.07 , a=1.0 in. 


<r/c. 
0 0 

Mi=M 4 

-*f*o * /<F o 

^=21*4 

-t /* 0 <r/<r 0 

1.716 .146 

1.852 0.0 

1.938 -.103 

.277 .0239 

.299 0.0 

.313 -.0169 

.142 .0128 

.153 0.0 

.161 -.00908 

.0 9 84 . 00 940 

.107 0.0 

.112 -.00674 

.0776 . 007 9 8 

.0841 0.0 

.0 8 83 -.0 0 577 

.0657 .00729 

.0714 0.0 

.0750 -.00532 

.0579 .00690 

.0631 0.0 

.0663 -.00508 

.0522 .00659 

.057 0 0.0 

.0600 -.004 90 

.0475 .00625 

.0520 0.0 

.0549 -.00470 

.0437 .00585 

.04 80 0.0 

.0506 -.00445 

.0406 . 0 0 53 8 

.0446 0.0 

.0472 -.00415 

.0381 .004 87 

.0420 0.0 

.0445 -.00380 

.0361 .00432 

.0399 0.0 

.0423 -.00341 

.0344 .00374 

.0380 0.0 

.0404 -.00300 

.0328 .00315 

.0364 0.0 

.0387 -.00257 

.0312 .00255 

.0347 0.0 

.0370 -.00213 

.0297 .00194 

.0330 0.0 

.0352 -.00167 

.0281 .00135 

.0313 0.0 

.0334 -.00121 

.0264 .00076 

.0295 0.0 

.0315 -.00075 

.0247 .00019 

.0277 0.0 

.02 96 -.00030 

.0230 -.00035 

.0258 0.0 

.0276 .00013 

.0212 -.00087 

.0239 0.0 

.02 56 . 000 56 

.0195 —00136 

.0219 0.0 

.0235 .00096 

.0176 -.00181 

.0199 0.0 

.0214 .00133 

.0140 -.00260 

.0158 0.0 

.0170 .00200 

.0122 -.002 93 

.0137 0.0 

.01-48 .0022 9 

.0103 -.00322 

.0117 0.0 

.0126 .00254 

.00657 -.00366 

.00744 0.0 

.00 803 . 002 92 

.0 0 2 82 -.0 0 3 91 

.00320 0.0 

.00345 .00314 

.00094 -.00396 

.00107 0.0 

.00115 .00318 
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Table 25. The effect of the material properties on the 

normalized stress intensity factors (k.j/k 0 , k 2 /k Q ) 
and the strain energy release rate, G/G 0 , with 
different thickness adherend3. 

(iya)= 1.0 , (iya)s0.5 » (hg/aJsO^ , 
ko-ToVC » G o =<r o a/E 2 » a=1 *° ±a - 




Hi=M 4 

|# 1 =2M 4 

k i/k 0 

-.00159 

-.00277 

-.00356 

k 2 /k 0 

-.0123 

-.0133 

-.0140 

G/Go 

• 399E-3 

.467E-3 

.517E-3 


Table 26. The effect of the material properties on the 

normalized stress intensity factors (k.j/k 0 , k 2 /k Q ) 
and the strain energy release rate, G/G Q , with 
different thickness adherends. 

(iya)=0.5 , (iya)=1.0 , (hg/aJsO.O? , 
k o=*o'F> G o =<r o a/E 2 > 3=1 *° in * 



P^Mij/2 

Pl=M4 

|» 1 32P2! 

k 1 /k o 

.00413 

.00277 

.00133 

k 2 /k o 

-.0122 

-.0133 

-.0142 

G/Gq 

.438E-3 

.48TE-3 

.537E-3 
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Table 27. The effect of the material properties on the 

normalized stress intensity factors (k.|/k 0 , k 2 /k 0 ) 
and the strain energy release rate, G/G QI with 
equal thickness adherends. 

OyaMiyaJsO.S , (h2/a)=0.07 , V«T 0 \/a 
Go =<r o a/E 2 • ^ ,0 in * 



P^Pj/2 

Mi=M 4 

M 1 =2P 1> 

k 1 /k o 

.00158 

.0 

-.00129 

k 2 /k c 

-.0106 

-.0118 

-.0150 

G/G 0 

• 303E-3 

•367E-3 

.5 94E-3 


Table 2 8. The effect of the material properties on the 

normalized stress intensity factors (k.|/k 0 , k 2 /k Q ) 
and the strain energy release rate, G/G 0 , with 
equal thickness adherends. 

(h-j/aJsOty/aJsl.O , (h2/a)=0.07 , k 0 =* 0 \fa 

G 0 =<ypa/E 2 , a=1.0 in. 



P 1 = Pl}/2 

Pl=P 4 


V k o 

.00113 

.0 

-.0008 

k 2 /ko 

-.0133 

-.0143 

-.0150 

g/Qq 

.470E-3 

.540E-3 

.594E-3 
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FIGURES 
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8 ' 





1: h«/a= 1.00 
2: h, /a=0. 50 
3: h^/asO.25 

Figure 4. The adhesive shear stress for different (h.j/a 
ratios with equal thickness adherends. 

(h ? /a=0.0025 , h.. = h,. , a=1.0 in. 




0.002 0.004 


0-006 0.008 

I 


Figure 5. The nornalized stress intensity fac 
strain energy release rate versus ( 
with equal thickness adherends. 

(.h^/a)=(h ii /a)=0.25 , a=1.0 in. 
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Figure 



h,/a 


6. The normalized stress intensity factor and the 
strain energy release rate versus (h-|/a) ratio 
with equal thickness adherends. 


(hp/a) =0.0025 , a=1.0 in. 



Figure 7. 


Figure 8. 



Comparison of the strain energy release rate 
calculated in this study (G) with the plate solutions 
(Gp)for the specific geometry discussed in Tbl. 9 . 



The effect of (h^/a) ratio on the normalized stress 
intensity factors with different thickness adherends. 
(h2j/a)=0.5 , (hg/aJrO.I , a=1 .0 in. 


93 




Figure 9- pie effect of (hg/a) ratio on the normalised stress 
intensity factors and the strain energy release rate 
with different thickness adherends. 

(h^/a)=1 .0 , (h^/aJsO.5 , a=1.0 in. 


S4 









0.5 1-0 2.0 

Figure 11. Effect of material properties on normalized stress 

intensity factors with different thickness adherends. 
(h^/a)=1.0 , (h^/aJrO.S , (hg/aJzQ.OE , a=1.0 in. 



0.5 10 2.0 fX,/^ 

Figure 12. Same as Figure 11 with 

( h./a) =0. 5 , ( hh/a) =1 . 0 , (h 2 /a)=0.07 , a=1.0 in. 
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APPENDIX I 


c = M 2 / 1^ i d = N 3 / Ri{ 

3 l ^ = (1^2 “ K^c) / (c “ 1 ) 1 

= (K^ C + 1) / (O - 1) , 

a 3 = ^ k 2 +c)/(°“D » 

= ( 3-2 ^2 ^ a 3^ ^ ) » 

b 1 = (k 3 - Kj,d) / (d - 1) , 

b 2 = (K^d + 1) / (d - 1) 
b 3 = (k 3 + d) / (d - 1) 

bjj = ^ b 2 h 3 + ^3^4^ » 

where P i , s, h^s, K^s are defined in Figure 1 
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APPENDIX II 


fi 

= e- 2<Xh 2 


fi* 

= e- 2Cth 3 


r 

= e- 2<X(h 1 

+ hg) 

*f * 

= e - 2<t(h 4 

+ h 3 ) 

* = 

: e- 2flth 1 


ii 

* 

: e - 2cth 4 



G^flC) = 2«Yaj, + 2CC(h 2 $ - fa^) , 

G 2 (CC) = a 2 Y + £ + + ^a 2 h 1 h 2 | - a 3 , 

G 3 (a) = - a 2 Y - fi + 20^? , 

Gjj(flt) = - 4a 2 h.|h 2 ,£ - 3.2 ~ 5 - a.j£ + a^Y > 

Gg(«) = 2a ai} + 2a ( h 2 5 - h^) , 

Gg(«) = 2a + $ + a 2 t 

g 7 («) = G^aJGgta) - G 3 (a)G 4 (o) , 

Gg(o) = G 3 (a)G 5 (a) - G 2 (a)G 6 (a) , 
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G 9 «X) = - 2«b 1< - 2CC(h 3 f* - h^ 8 ) , 

G 10 (a) = b 2 + $* + i4« 2 h 3 h^* + b^ 8 - b 3 Y 8 , 
G n (a) = - b 2 - | 8 - 2«th lj * 8 , 

G 12 («) = - b^* + b 3 - iiOC 2 h 3 h J4 f * - fi* - b 2 Y* , 
G 13 («) = - 2Ct*f *bjj - 2a(h 3 ^ § - h 4 f*) , 

G u (a) = b 2 Y* - 2«h 4 f 8 + fi* , 

G 15 (a) = G 11} (ct)G 9 (a) - G^CcOG^Ca) , 

G 1g (a) = G^CcOG^Ca) - G u (a)G 10 (a) , 

where h^'s, a^'s, b^’s are defined m Appendix I . 
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APPENDIX III 


P^cc) = G^a) , 

P 2 («) = c 2 (a) , 

P 3 (CC) = -a 2 Y - £ - zahfi , 

P 4 «n) = g 4 (o) , 

P 5 (ct) = G 5 (a) , 

P 6 (a) = 2«h^ - I - ^ , 

p ? (a) = p^cOPgCa) - P 3 (a)P 4 (a) , 

P 8 (a) = P 3 («)P 5 («) - P 2 («)P 6 (a) , 

P 9 (a) = G g (a) f 

p 10 «x) = G 10 (a) , 

P^ (00 = - b 2 - |* + 2<xh^ , 

P 12 (a) = G 12 (a) , 

p 13 («) = g 13 «x) , 

P-lij(ct) = - b 2 Y* - 2<xh iJ $* - £* , 

P i5 (a) = P 11} (a)P 9 (a) - P 12 (cc)P n (a) 
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P l6 (ct) = p^tflOP^tct) - P 12j (a)P 10 (a) , 


where (h^s, a^s, b^s) and (fi , £*, Y, Y*, 5 , ?*, G i , 3) are 
defined in Figure 1 and Appendix II, respectively. 
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APPENDIX IV 


Solution of the Singular Integral Equation 


In [15] a quadrature formula of closed type is derived for the 
principal part of the singular integral equation of the form ; 


a #(t) + — 
-I 


♦ (t) 


t-x 


I 

-/■ 


dt + | ♦(t) k(x,t) dt = g(x) , 

-I 

-1 < x < +1 


The solution will be sought in the form 


♦ (t) = *(t) W(t) , 


where ^(t) is a bounded function, and W(t) is given by 


W(t) = (1-t) a (1+t) 


* 


1 

« = log 

ZVi 


r a - ib 


a + ib J 


+ N 


1 

log 

2ir i 


a - ib 


L a + ib J 


+ M . 


105 



N and M are arbitrary integers determined from the physics of the 
problem. The index of the singular integral equation is defined by 

K S - (« + fi) = - (N + M) . 

Defining, 

, b 

k (x,t) = + k(x,t) , 

ir(t-x) 

the approximate solution to the singular integral equation is 
determined from 

■rv-i 

(1 + O) H 1 k*(x k ,t.,) \ k *( x k» fc i^ 2 ) 

Ua, 

+ <1 + *) Hq k # (x k ,t n ) *(t Q ) = g(x k ) , k=1 . ,-1 

and 

*v-t 

(1 + «) H 1 *(t 1 ) + E H i ^ (t i } + (1 + *) H n 4»(t n ) = C , 
where the last equation is the approximation to the extra condition, 
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(when K=1) 


l 

/♦( t) dt = C . 

-I 

Here t lf and H i are defined as follows , 


(1 - t|) 1+ * } (t ± ) = 0 , t, > t 2 >....> t n 

p £zj- a ’ “ 1 "^ ) (x k ) = 0 , k=1 ,2 . , (n-1 ) 


and 


T(n+a) I*(n+£) 

H± " (n-1) [r(n)P^«^> (t ± )] 2 

where P^® ^ (x) is the Jacobi polynomial of degree "n" . 

In the analysis of the report the singular integral equations 
are of the first kind and CC= £= - 1/2 giving K = 1 . In this case 
Jacobi polynomials reduce to the Chebyshev polyncnials and n 
unknowns 4>(t^) can be determined from 
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1 «-» 

— k*(x k ,t.,) «I»(t.j) + Y. k*(x k ,t.j 
2 ‘ 


a.«a 


+ -k*(x t ) *(t n ) = 

2 *0 


and, 


1 1 
— *(t,) +1 ♦ (tj.) + -«Kt Q ) 
2 L-a. 2 


where 


= Cos 


i-1 


L n-1 


, i= 


x^ = Cos 


2k-1 


L 2n-2 


» k: 


) *(t ± ) 


g(x k ) , 


n-1 

C , 

■w 


!}•••• j n 


!>•••• ?( n- 1 ) 
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